Centre de Recherches Mathematiques 
CRM Proceedings and Lecture Notes 

2000]Primary 34C07, 34C08, 34M35; Secondary 34C10, 34M15, 34M50, 14Q20, 
32S65, 13E05 

Quantitative theory of ordinary differential equations 
and tangential Hilbert IG**^ problem 

Sergei Yakovenko 

Abstract. These highly informal lecture notes aim at introducing and ex- 
plaining several closely related problems on zeros of analytic functions defined 
by ordinary differential equations and systems of such equations. The main 
incentive for this study was its potential application to the tangential Hilbert 
16th problem on zeros of complete Abelian integrals. 

The exposition consists mostly of examples illustrating various phenom- 
ena related to this problem. Sometimes these examples give an insight concern- 
ing the proofs, though the complete exposition of the latter is mostly relegated 
to separate expositions. 



Lecture 1. Hilbert 16th problem: Limit cycles, cyclicity, Abelian 
integrals 

In the first lecture we discuss several possible relaxed formulations of the Hilbert 
16th problem on limit cycles of vector fields and related questions from analytic 
functions theory. 

1. Zeros of analytic functions 

The introductory section presents several possible formulations of the question 
about the number of zeros of a function of one variable. All functions below are 
either real or complex analytic in their domains, eventually exhibiting singularities 
on their boundaries. We would like to stress that only isolated zeros of such func- 
tions are counted, so that by definition a function identically vanishing on an open 
set, has no isolated zeros there. 

Exposition goes mostly by examples that are separated from each other by the 
symbol -4. A few demonstrations terminate by the usual symbol □. 

1.1. Nonaccumulation and individual finiteness. A function f{t) real 
analytic on a finite open interval (a, 6) C M may have an infinite number of isolated 
zeros on this interval only if they accumulate to the boundary points a, b of the 
latter. Thus the finiteness problem of decision whether or not the given function 
/ has only finitely many zeros in its domain, is reduced to studying the boundary 
behavior of /. 

In particular, if / is analytic also at the boundary points a, b, then accumulation 
of infinitely many zeros to these points is impossible and hence / has only finitely 
many roots on the interval (a, b). However, this strong condition of analyticity can 
be relaxed very considerably. 
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Example 1.1. Assume that f{t) defined on (0,1] admits a nontrivial asymp- 
totic expansion of the form 

oo 

f(t) r^Y' Ckf" asi^O+, 

h (1-1) 

Cfe, Tfe e M, ro < ri < • • • < rfc < • • • 

i.e., the difference between / and a partial sum decreases faster than the next 
remaining term (the nontriviaUty means that not all Ck arc zeros). Then in a 
sufficiently small semi-neighborhood of i = the function / has the same sign as 
the first nonzero coefficient Ck (since the functions f are nonvanishing on M_(_), and 
hence roots of / cannot accumulate to implying that / has only finitely many 
zeros on (0, 1]. M 

This example can be generalized for functions admitting asymptotic expansion 
in any system of mutually comparable functions /i,/2,..., when fk+i = o{fk), 
provided that all of them keep constant sign near the boundary point (s). 

Example 1.2. Any function on (0,a], a > 0, representable as a finite sum of 
the form 

hkr{t) f In''-^ t, r G R, A; G N, 

fe,r 

with the functions hkr{t) real analytic on [0, a] (i.e., including the boundary t = 0), 
cannot have infinitely many roots accumulating to t = 0. Indeed, the above sum 
expands in the monomials In'^ r G M, A; G Z, that are naturally lexicograph- 
ically ordered by their growth rates as i — > 0+. Clearly, a function admitting an 
asymptotic expansion involving terms of such form, also possesses the finiteness 
property. 

In this example we do not exclude the cases when the expansion is trivial. 

However, the convergence assumption on h^r implies than in such cases the function 
is identically zero and has no isolated roots at all. 

1.2. Parametric families of analytic functions, localization and cycl- 
icity. Consider a function / = f{t, A) = f{t, Ai, . . . , A„) real analytic in an open 

domain U of the space K x K" (this means that / can be expanded in a converging 
Taylor series centered around any point in [/). The function / can be considered 
as an analytic family of functions fx = f{-,X) defined in variable domains U\ = 
Un(Rx {A}) C M depending on A. 

Our nearest goal is to formulate the parametric finiteness property and establish 
simple sufficient conditions for it, similar to the non-parametric case. 

Definition 1. Let A be a point set in K or C. Everywhere below we denote 
by #j4 the number of isolated points (finite or not) of A. 

Definition 2. We say that the analytic family / = {fx}: U ^ R possesses 
the uniform finiteness property, if the number of isolated zeros of all functions fx 
in their respective domains Ux is uniformly bounded by a constant independent of 
A. 

In the same way as in the non-parametric context, it would be desirable to 
derive the uniform finiteness from some local properties of the family /. 
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Definition 3. The cyclicity of the family / at a point (t*, A*) e M x K" is the 
upper limit (finite or not) 

^/(t,,A,) =limsup sup #{teUx:\t-U\<s, fxit) =0}. (1.2) 

e^0+ ||A-A.||<£ 

The term comes from the bifurcation theory (see below). 

We stress again that the notation ■■} above means the number of isolated 
roots of fx; the upper limit may be infinite even if all these numbers are finite. 

If fx, ^ 0, then cyclicity of any family containing this function can be majorized 
in terms of only the function itself. 

Example 1.3. Suppose that / is analytic at an interior point (0,0) e U, and 
/o ^ 0. Then the cyclicity of / at the origin is finite. Moreover, if t = is an 

isolated root of multiplicity /z, that is, f{t) = cf^ -\ , 0, then the cychcity is 

no greater than /x, that is, A//(0,0) ^ /x. 

This follows from the fact that the /ith derivative of / is nonvanishing at the 
origin and hence at all sufficiently close points. But a function whose derivative 
/(")(*) has a constant sign on an interval, cannot have more than /z isolated zeros 
(even counted with multiplicities), as follows from the iterated RoUe theorem. 

This example can be easily generalized for the case of complex analytic functions 
defined in a domain U c C x C". In the above assumption /q**^ (0) ^ 0, the cyclicity 
of the family will be exactly ji. To prove this, one can choose a small circle around 
the origin and apply the Rouche theorem to it. Another possibility would be to use 
the Weierstrass preparation theorem. < 

Example 1.4. If dim A = 1, then any function analytic in (f. A) near t = tit, 
A* = can be expanded as f{t, A) = /o(t) + \fi{t) + \^ f2{t) + ■■■■ If f{t, A) ^ 0, 
then for some finite k necessarily fk{t) ^ 0, and after division by A*^ the question 
about isolated zeros of / can be reduced to the situation when f{t, 0) ^ 0, discussed 
earlier. In this case A/'/(f*,0) is no greater than the multiplicity of fk at f*, where 
fk is the first nonzero term in the expansion. < 

This example illustrates an absolutely general fact about analytic functions (no 
matter, real or complex): the cyclicity M}{t,\) takes finite values at all interior 
points of the domain of analyticity U. 

Theorem 1. /// is analytic at a point {t,\) e R x M", then Mf{t,\) < +oo. 

This assertion can be derived from general finiteness properties of analytic sets, 
see [Loj91]. In the present form the theorem was formulated in connection with 
bifurcations of limit cycles, see [Rou98] . 

The proof of finiteness of A// in the general multiparametric case requires some 
analytic techniques. One possibility — assuming for simplicity the point to be at 
the origin (0, 0) — is to consider the expansion f{t, A) = X^^^o ak{\)t'^ and the ideal 
in the ring EH of analytic germs at (K" , 0) , generated by the coefficients ak (A) . 
This ideal is called the Bautin ideal [Rou98]. Since the ring D\ is Noetherian, 
the Bautin ideal is in fact generated by a finite number of the coefficients ak- If 
the germs ai,...,a,y generate the Bautin ideal for some finite i^, then one can 
show that the function f{t, A) can be represented as Ylk=o <^k{^)t''hk{t, A) with hk 
analytic in t, A and hk{0, 0) ^ 0, see [Rou98]. Prom this it is already easy to derive 
that A// (0,0) ^ u < oo. Almost no modification is required to cover the complex 
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analytic case as well: the only difference is that the usual Rolle cannot be applied 
to holomorphic non-real functions. The alternative is to use the complex Rolle 
theorem from [KY96] . An example of such use in a situation very similar to that 
discussed above, can be found in [YakOO]. 

1.3. Finite cyclicity and uniform finiteness. The main reason for intro- 
ducing the notion of cyclicity is the following very simple but basic theorem. In 
applications to bifurcation of limit cycles of analytic vector fields, it was observed 
by R. Roussarie [Rou88, Rou89, Rou98]. 

Let / : [/ ^ M as before be an analytic family of functions and A// : U — > 

N U {+00} its cyclicity function. 

Theorem 2. // the closure U is a compact subset oJRx M", and cyclicity Aff 
is finite everywhere on U, then the family f = {f\} admits a uniform upper hound 
on the number of isolated roots: 

sup #{t e C/a : fx{t) = 0} < +00. (1.3) 

A 

Proof. By definition of the counting function, any point (t,. A,) e U can be 
covered by a sufficiently small cube {|t — t*| < e, 1|A — A*]] < e} such that the 
number of isolated roots of / in this cube is no greater than the number A^j(f», A*) 
finite by the assumptions of the theorem. It remains only to choose a finite subcover 
of the compact U and add together the corresponding cyclicities. □ 

Theorem 1 asserts that the cyclicity is automatically finite at all interior points 
of the domain i7, so it is again the boundary behavior of a given parametric family 
that determines whether this family possesses the uniform finiteness or not. 

1.4. Terminology: individual vs. existential finiteness problems. The 

arguments proving finiteness of the number of isolated zeros of an individual an- 
alytic function f{t) and the arguments establishing uniform boundedness of this 
number for parametric families of such functions, are both of the same purely ex- 
istential natiu'e: neither of them gives any way to com,pute or even estim,a,te from, 
above these numbers. However, the parametric claim is definitely stronger than the 
assertion concerning individual functions of this family (one can easily construct 
families in which the number of isolated zeros is always finite but not uniformly 
bounded). 

We shall repeatedly encounter the problem on bounding the number of isolated 

zeros for various classes of analytic functions and finite parameter families of such 
functions, mostly defined by ordinary differential equations with polynomial right 
hand sides. For each class one can pose several finiteness- type problems in the 
increasing order of strength (the gap in the enumeration will be explained below) . 

1. Individual finiteness problem. Prove that each function from the family 
possesses only finitely many isolated zeros. 

2. Existential {^iniform) finiteness problem. Prove that the number of isolated 
roots is uniformly bounded over all functions from this family. 

4. Constructive finiteness problem. Find an explicit upper bound for the num- 
ber of isolated roots or at least find an algorithm for computing this bound. 
The adjective "existential" stresses the difference between the last two types of 
problems, whereas the adjective "uniform" will be used to underscore the difference 
between the first two assertions. 
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Example 1.5. Consider the class of parametric families of functions real an- 
alytic on [0, 1] depending analytically on parameters from [0, 1]" C R"+^ (in both 
cases including the boundary). As was just explained, the individual finiteness for 
functions from this class follows from the uniqueness theorem for analj^ic func- 
tions while the existential (uniform) finiteness theorem (Theorem 2) follows from 
Theorem 1. 

Example 1.6. Polynomials in one variable of degree ^ d for any finite d form 
a finite-parameter family. Constructive finiteness theorem for this "toy" class of 
functions is known as the Principal theorem of Algebra. Less trivial examples can 
be found in §8. M 

1.5. Constructive finiteness. In practice the existential finiteness is always 
derived from finite cyclicity using Theorem 2. Both steps (finite cyclicity and its 

globalization) use arguments of existential nature. However, at least theoretically it 
may happen that the function A// is explicitly bounded at all points of the closure 
U. This still would not allow to compute explicitly the global uniform bound 
on the number of isolated zeros of /, but such local constructive finiteness would 
be clearly a much stronger assertion concerning the family /. The corresponding 
finiteness problem occupies an intermediate place between existential and (global) 
constructive finitc^ness problems. 

3. Constructive finite cyclicity problem. Find an explicit majorant for the 
counting function A / . 

Yet it should be remarked that in order to discuss the constructive finiteness prob- 
lems (global or local), the family should be defined by some algebraic data, other- 
wise computability does not make sense. We postpone discussion of these question 
until §8. 

Example 1.7. Let Fi{t), . . . , Fn{t) be analytic functions (real or complex), sat- 
isfying together a system of polynomial ordinary differential equations 
dP 

Y. '^ikat'F'^, F« = Ff^...j;«", i = l,...,n, (1.4) 

k+\a\^d 

assuming the degree d and all coefficients Cika (real or complex) known. Then 
for any finite m G N the polynomial combinations of these functions constitute a 
finite-parametric family: the coefficients A = {Xka : |q;| ^ to} of the polynomial 
combinations /(t. A) = J2k+\a\^m, ^i"^ ^''■^"i^) being the parameters. The domain 
of this family is clearly the Cartesian product of the linear space of the appropriate 
dimension and the common domain U of the functions Fi. Though solutions of 
polynomial systems may blow up in finite time, we will assume for simplicity that 
all Fi are analytic on the compact closure U of U. 

In these assumptions both individual and existential finiteness are guaranteed 
by Theorem 2 which, as usual, does not provide any bound on the number of 
zeros. Yet using algebraicity of the system defining these functions, one can derive 
explicit bounds. In [Gab99] A. Gabrielov found an upper bound for the maximal 
order of zero {multiplicity of a root) that a nontrivial (not identically equal to zero) 
function from the family may have, in terms of n,d and m. The answer is 
given by an expression polynomial in the degrees d, m and exponential in n (the 
dimension of the system). Y. Yomdin in [Yom98] derived from this multiplicity 
bound the corresponding cyclicity hound, using arguments involving Bautin ideal 



6 



SERGEI YAKOVENKO 



and its generators (a subtle generalization of the approach outlined in discussion 
of Theorem 1). This combination of the two results gives a complete solution for 
the constructive finite cyclicity problem for polynomial combinations of functions 
defined by systems of polynomial ordinary differential equations. M 

The corresponding global problem was settled by D. Novikov and the author 

[NY97, NY99b]. Unlike the bound for cyclicity, the bound for the total number 
of zeros depends on the size of the domain p = maxtgt/ \t\ where they are counted 
and, rather naturally, on the magnitude of the coefficients R = maxj \cika\ of 
the system (1.4). The answer is polynomial in p, R but as a hmction of the other 
(integer) parameters n, to, d is a tower of height 4. The detailed explanation of this 
result is postponed until §11. 

1.6. Paradigm. All the above results on elementary properties of analytic 

functions and their parametric families, were mentioned in order to set out a para- 
digm that will be used when discussing much more complicated parallel problems 
on limit cycles of analytic and polynomial vector fields. In order to make the sim- 
ilarity transparent, all finiteness theorems can be organized in the form of a table 
with finiteness types labeling the rows and corresponding theories represented by 
different columns. The first column describes theory of analytic functions as in 
Example 1.5 above. Cells of this column list the key theorems ensuring each type 
of finiteness types. 

The following section §2 describes some known results fitting the second column 

of Table 1: the strongest form of the Hilbert 16th problem should occupy the 
place where the question marks appear. The parallel non-constructive results for 
Abelian integrals will be discussed in §4: their complexity level can be described 
as intermediate between relatively elementary theorems on analj^tic functions and 
transcendentally difficult counterparts on limit cycles. 

Remark. In general, the results in the lower part of the table are more dif- 
ficult than those from the upper part. However, this is not true when comparing 
existential finiteness with constructive finite cyclicity: while the former deals with 
boundary behavior, the latter addresses the issues of completely different nature 
that can be (and indeed sometimes are) simpler. In particular, this is the case 
with Abelian integrals, where the bound for multiplicity can be be relatively easily 
obtained, see [MarQl] and §6.5 below. 

2. Limit cycles of planar vector fields 

This section briefly surveys a few known general finiteness results for limit 

cycles of vector fields. 

2.1. Beisic facts. A polynomial differential equation 

x = P{x,y), y^Q{x,y), P,QeR[x,y] (2.1) 

defines a vector field on the real plane M-^, spanning the distribution (line field) 
with singularities given by null spaces {ui = 0} of the polynomial Pfaffian form 
ft = Q dx — Pdy. If P and Q have a nontrivial greatest common divisor R, then 
this distribution on the Zariski open set {R ^ 0} coincides with that given by the 
form R~^il, and hence extends analytically to all but finitely many points of the 
curve {R = 0}. For the same reasons the distribution extends onto the infinite line 



QUANTITATIVE THEORY OF ODE 



7 



Finiteness 
type 


Category of objects 


Functions 
analytic 
on a compact 


Limit cycles 
of polynomial 
vector fields 


Abelian 
integrals 


Individual 
finiteness 


uniqueness 
theorem 


Ilyashenko (1991)- 
Ecalle (1992) 


easy exercise 
see §4.1 


Existential 
finiteness 


finite cyclicity 

Theorems 1 
and 2, sec §1.3 


in progress for 

quadratic case (121 
polycyclcs [DRR94]) 


Varchenko- 

Khovanskii 
(1985) 


Constructive 
finite 
cyclicity 


Gabriclov- 
Yomdin (1998) 


known for fc-gcneric 
elementary and few 
other polycycles 


Mardcsic 
[Mar91] 
(multiplicity) 


Constructive 
finiteness 


Novikov- 
Yakovenko 
(1999) 


??? 

(nothing in view) 


see Table 2 



Table 1 . Various flavors of the Hilbert 16th problem 



MP^ on the projective plane MP^, having at most finitely many isolated singulari- 
ties on the latter after cancellation of an eventual common factor. In other words, 
when dealing with an individual differential equation (2.1), one can assume that the 
corresponding distribution has only isolated singularities. It is sometimes conve- 
nient to talk about foliations with isolated singularities defined by the distribution 
{n = 0} on MP2. 

A limit cycle 7 by definition is an isolated periodic trajectory of (2.1). Limit 
cycles of a vector field can be tracked by their intersections with analytic arcs 
transversal to the field. More precisely, near any intersection p = 7 n cr of a periodic 
orbit 7 with a transversal arc a, the Poincare first return map A^-, the transport 
along integral curves until the next intersection with the transversal, is defined. It 
has p as a fixed point, isolated if 7 is a limit cycle. Choosing an analytic chart t 
on a allows to describe p as a root of the displacement function f{t) = Acr{t) — t. 
All nearby limit cycles must intersect a by points that are isolated roots of the 
displacement /. 

From the general theorem on analytic dependence of solutions of differential 
equations on initial conditions and parameters, the displacement function / is an- 
alytic near p. As a corollary, we conclude that an infinite number of limit cycles 
(corresponding to isolated roots of the displacement function) cannot accumulate 
to a periodic orbit of the field. In a similar way, if the vector field depends ana- 
lytically on parameters A and exhibits a periodic orbit 7 (isolated or not) for one 
value A* of the parameters, then there exists a finite upper bound for the number 
of limit cycles in a small annulus around 7, uniform over all values of A sufficiently 
close to A*. 

In general, one transversal cannot "serve" all limit cycles, and there is no nat- 
ural way to define the maximal domain of the Poincare return map unambiguously 
and globally. However, the discussion in §1 suggests that it is the boundary behavior 
that is important for counting isolated zeros of the displacement. 
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2.2. Polycycles and limit periodic sets. Instead of trying to reduce for- 
mally the global investigation of limit cycles to that of return maps for one or several 
transversal arcs, it is better to study the cycles (compact leaves of the foliation) 
themselves. 

For the non-parametric case, the analogue of a boundary point of the domain 
of the Poincare return map is a polycycle, an invariant set consisting of one or 
more singular points of the vector field, and a number of bi-infinite trajectories 
connecting them in a cyclic order (repetitions of singular points allowed). More 
accurately, one can show that the only sets that can appear as Hausdorff limits of 
periodic orbits of a vector field having only isolated singular points, are polycycles. 

In the parametric case one cannot assume anymore that singularities of the 
foliation are isolated (the polynomials P and Q may have common factors for some 
values of the parameters and be mutually prime for the rest). Still one can show 
[FP86] that the Hausdorff limit of a family of limit cycles occurring for converging 
values of the parameters, must be a limit periodic set (also known as graphic), 
the object differing from a polycycle in only one instance, namely, it may contain 
analytic arcs of non-isolated singularities. 

2.3. Individual finiteness theorem and finite cyclicity. Following the 
paradigm set out in §1.3, one can easily formulate counterparts of the individual 
and existential finiteness problems for limit cycles for vector fields. The individ- 
ual finiteness theorem (known also as the Dulac conjecture after being believed for 
some 60 years to be the Dulac theorem) asserts that any polynomial vector field can 
have only finitely many limit cycles. This is an easy corollary to the highly nontriv- 
ial nonaccum,ulation theorem, asserting that limit cycles of an analytic vector field 
cannot accumulate to any polycycle. The latter assertion is the most spectacular 
and the most general fact established so far in connection with the Hilbert 16th 
problem, that was independently proved by Ilyashenko [Ily91] and Ecalle [Eca92] 
by totally different methods. This is a typical example of an assertion on boundary 
behavior of the Poincare map. 

Remark. It would be appropriate to notice here that Dulac reduced the in- 
dividual finiteness problem to the situation discussed in Example 1.1. However, 
he did not notice that the expansion (1.1) (only involving logarithms) may well 
be trivial, hence not allowing for so easy treatment. It took both Ilyashenko and 
Ecalle hundreds of pages to prove that even in this apparently highly degenerate 
case the nonaccumulation still holds. 

The existential finiteness problem for polynomial vector fields arises very natu- 
rally, since coefficients of the polynomials P, Q can be treated as parameters. These 
parameters can be considered as ranging over the compact projective space of a siiit- 
able dimension, since simultaneous multiplication of both P and Q by a nonzero 
constant does not alter the phase portrait. Existence of a uniform upper bound for 
the number of limit cycles would mean that for any degree d the number of limit 
cycles for a polynomial vector field of degree ^ d is bounded by some number 'H{d) 
depending only on d. According to the Roussarie localization theorem (an analog 
of Theorem 2 for limit cycles, see [Rou88]), it would be sufficient to prove that 
any limit periodic set has finite cyclicity within the universal family of polynomial 
vector fields of degree < d. 
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However, even for the simplest nontrivial case of quadratic {d = 2) vector fields, 
this is not proved. Dumorticr, Roussaric and Rousseau composed in [DRR94] a 
list of 121 graphics occurring for quadratic vector fields, and reduced the existential 
finiteness problem to proving that all these graphics have finite cyclicity. First 
several cases were studied in the same paper by applying relatively standard tools 
of bifurcation theory. Since then many other considerably more delicate cases 
were investigated and their finite cyclicity proved, so that one can hope that the 
existential finiteness will be ultimately established in the qiiadratic case. At the 
same time it is clear that such type of case study is even theoretically impossible 
for higher degree cases. 

On this background the perspective of computing or even explicitly m,ajorizing 
the ^^Hilbert number" which is the strongest (constructive) form of the initial 

Hubert's question "on the number and position of limit cycles" [HilOO], looks very 
remote. 

2.4. Digression: constructive solutions of localized problems. As was 

already noted above, though finiteness of cyclicity for all graphics is not proved, 
one can independently work towards obtaining constructive bounds for cyclicity in 
cases when its finiteness is already known. This is possible since cyclicity is very 
strongly depending on the types of singular points lying on the polycycle (resp., 
the limit periodic set). 

For example, if all singular points on the polycycle are elementary (having 
non-nilpotent linear parts), then cyclicity of generic r-parametric families can be 
estimated in terms of the number r of independent parameters. The corresponding 
algorithm was suggested by Ilyashenko and the author [IY95b] and improved to a 
very concise and explicit answer by Kaloshin [IK99] . Some other cases of polycycles 
(usually carrying one or at most two singular points) also have known cyclicity, see 
[KS95] and references therein for a synopsis. These results are partly reflected in 
Table 1. 

Yet if the problem is formulated for polynomial vector fields and the bounds 
are required to be given in terms of the degree, the problem immediately becomes 
transcendentally difficult. Almost nothing is known, for example, on the maximal 
multiplicity of a (nonsingular) limit cycle of a polynomial vector field. Cyclicity of 
singular points is a considerably more "algebraic" problem. For example, when the 
linear part of the singularity is a nondegenerate rotation, a polynomial algorithm 
can be suggested which stops at a step number A'' if and only if the cyclicity of 
a singular point is N. However, the running time of this algorithm is absolutely 
unknown (see the discussion of a similar problem in §3.3). Recently M. Briskin, 
J. -P. Frangoise and Y. Yomdin treated in details one problem on Abel equations 
that may give a clue to understanding some of the phenomena, see e.g., [BFY98a, 
BFY98b]. 

The following section describes one very important case when the local (with 
respect to parameters) constructive finiteness can be reduced to an algebraic con- 
text. 

3. Abelian integrals: appearance and basic properties 

Consider a family of analytic vector fields X\ analytically depending on param- 
eters A e (M^jO) (we consider only a sufficiently small neighborhood of the origin 
in the parameter space). Assume that for A = the vector field Xq possesses a 
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limit cycle 70 of multiplicity < +00. As illustrated by Example 1.3, the cyclicity 

of 7o in this case is at most /i, being thus explicitly bounded. 

The problem becomes less trivial if the field Xq possesses an annulus filled by 
nonisolated periodic orbits, 70 being one of them. Vector fields with this property 
arc called integrablc or conservative. Such behavior, very unlikely for an arbitrary 
analytic family, becomes an event of finite codimension when polynomial vector 
fields are considered. 

To study bifurcations of limit cycles in perturbations of conservative systems, 
an approach described in Example 1.4 can be used. 

3.1. Perturbation of polynomial vector fields. Assume that the form 
fl = Qdx — Pdy defining the polynomial distribution {O = 0}, is exact: Cl = dH, 
where H = H{x,y) is a bivariate polynomial in two variables. The corresponding 
vector field is then a Hamiltonian one. 



with H being its Hamiltonian function or simply the Hamiltonian. 

Since dH vanishes on the null spaces of the distribution {f2 = 0}, each leaf of 
the corresponding foliation belongs to a level curve of H. In particular, a periodic 
orbit must be a compact oval of some level curve and hence all nearby close leaves 
must be also closed. Hence a Hamiltonian vector field cannot have limit cycles: 
the corresponding Poincare return map is identity and the displacement identically 
zero. 

Consider the following one-parameter family of polynomial distributions, per- 
turbing the Hamiltonian distribution: 

dH + suj = 0, s € {^^,0), LJ =pdx + qdy, p,q €R[x,y]. (3.2) 

The polynomial 1-form uj can be arbitrary. 

Consider an analytic segment a transversal to the oval 7(t*) C {H = t^}. This 
transversality implies that the Hamiltonian H restricted on <t, gives an analytic 
chart t on the latter. On the other hand, for all t sufficiently close to t^, one 
can unambiguously choose an oval ^{t) C {H = t} so that 7(i) tends to 7(i*) as 
t ^ tt:, e.g., in the sense of Hausdorff distance (this choice in a more broad context 
is discussed below, in §5). 

The following elementary computation gives the first variation of the displace- 
ment function f{t, e) = A„{t) — t with respect to the small parameter £ at £ = 0. 

Lemma 1 (Poincare-Pontryagin) . 



Proof. Consider the leaf (analytic curve) of the foliation {dH + euj = 0} 
passing through the point t on the transversal, and denote by ^{t, e) the (oriented) 
segment of this leaf between the the initial point and the next intersection with a. 
By definition of the chart t, the displacement function measured in the chart t, is 
the difference of H between the endpoints of 7(t, e), hence 




(3.3) 




dH 



QUANTITATIVE THEORY OF ODE 



11 



(the equality is exact). The form dH + eu vanishes on j{t, s), therefore the integral 

above is equal to the integral of — ew along "f{t,s) and, since j{t,s) converges 
uniformly to the dosed curve j{t) — ^{t, 0) C {H = t} as e ^ 0, we conclude that 

f{t, s) = —e (p u) + o(e) as £ ^ 0, 

with o(e) uniform and analytic in t and e. This yields the formula (3.3) for the 
derivative. □ 

Recall that an Abelian integral is the integral of a polynomial 1-form along 
any algebraic oval (i.e., the closed oval of a real algebraic curve). Clearly, Abelian 
integrals can be considered as functions of the parameters, the coefficients of the 
form and the algebraic equation defining the curve. However, in applications to the 
bifurcation theory, one coefBcient is clearly distinguished and plays the role of the 
argument, while others are treated as "true" parameters. Making a minor abuse 
of language, we shall always consider Abelian integrals as analytic functions in the 
following sense. 

Definition 4. A real Abelian integral corresponding to a Hamiltonian H and 

a polynomial 1-form lu = p{x, y) dx + g(x, y) dy, H,p,q G M.[x, y], is a multivalued 
function I{t) of the real variable t defined by integration of co over a real oval ^{t) 
of the algebraic curve {H{x, y) = t} C M?. 

The reason for multivaluedness is obvious, since there can be several ovals lying 
on the same level curve {H = t}. In §5 we describe the branching of the Abelian 
integral I{t) after analytic continuation into the complex domain. This complexifi- 
cation of Abelian integrals explains connections between different branches of these 
functions, see §5. 

3.2. Bifurcation of limit cycles from periodic orbits of Hamiltonian 
systems. Limit cycles appearing in the one-parametric system (3.2), correspond 
to isolated zeros of the displacement function f{t, s), real analytic in t and s, whose 
first Taylor term was just computed: 

fit,s) = Ioit)+sh{t)+s^h{t) + --- , Io{t) = 0, hit) = - <f to. (3.4) 

J-rit) 

This is exactly the situation treated in Example 1.4: if the first variation is non- 
trivial, Ii{t) ^ 0, then the number of isolated roots of f{t,e) on any interval of 
analyticity (i.e., not containing singularities at the endpoints) is no greater than 
the number of (necessary isolated) roots of the Abelian integral Ii defined by the 
polynomial H and the polynomial 1-form oj. Computing or majorizing the number 
of isolated zeros of Abelian integrals is the central theme of these notes: under the 
name Tangential Hilbert problem it is discussed below. 

The first variation method does not work if Ii{t) = 0. In this case higher 
variations /;s(')i ^ = 2,3,..., have to be computed and analyzed. 

The computation of higher variations is relatively simple. It follows from results 
of Ilyaslicnko [Ily69] and Gavrilov [Gav98] that for almost all Hamiltonians H (and 
certainly for Hamiltonians with isolated critical points, pairwise different critical 
values and transversal to infinity as defined below in §6), the condition /i = 
implies that there exist two polynomials G,F G R[x, y] such that 

w = G{x, y) dH + dF{x, y), hence du = dG A dH. (3.5) 
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(the inverse statement is obvious). 

Example 3.1 (see [Pra96]). Consider the Hamiltonian H{x,y) = ^{x'^ + y^). 

In the complex coordinates z = x -\- iy, z = x — iy wc have H = ■^zz. One can 
easily verify that the polynomial 1-form w = A{z, z) dz + B[z, z) dz has identically 
vanishing integral over the circles {H = t} if and only if the differential dui = 
{—Az + Bz) dz/\dz contains no monomial terms of the form (zz)'^ dzAdz. Any other 
monomial can obviously be represented in the form (3.5): z^z^ dz A dz = dG A dH 
withG = z'zV(«-i)- < 

Using the representation (3.5), the formula for the second variation of the 
displacement map can be expressed as an Abelian integral again. 

Lemma 2 ([Yak95, Pra96]). Assume that for the perturbation (3.2) the first 

variation vanishes identically, so that 



h{t) = j 



uj = 0, 

lit) 

and G is any polynomial satisfying the condition (3.5). Then the second variation 
of the displacement map is given by the following Abelian integral, 

h{t)=^ f{t,e) = -<f Guj. (3.6) 

This construction can be further iterated [Pra96] as long as necessary to ob- 
tain a variation that is not identically vanishing. If, on the other hand, all integrals 
obtained in this recurrent process, turn out to be identically zero, then the displace- 
ment function itself is identically zero and hence the perturbation (3.2) consists of 
integrable systems for all s. 

Remark. The exact conditions on the Hamiltonian H for the identity (3.5) to 
hold for any form whose integral is identically zero, involve connectedness of the 
complex afEne level curves H~^{t) for almost alH e C [Gav98]. Discussion of the 
multidimensional situation can be found in [BDOO]. As soon as these conditions 
fail (even because of coinciding critical values, for example), then computation of 
higher variations may involve integration of non-polynomial forms, see [Ili96]. 

3.3. Open problems: bifurcation of limit cycles from Hamiltonian 

polycycles and generalized Poincare center focus problem. It would be 
wrong to think that the above approach based on computing consecutive variations, 
completely reduces the question on limit cycles born by perturbation of Hamiltonian 
systems, to investigation of Abelian integrals. 

First, even assuming the simplest case Ii ^ 0, we cannot in general say anything 
about limit cycles born from critical level curves (corresponding to polycycles of the 
unperturbed Hamiltonian system). In the particular case of separatrix loops (level 
curves homeomorphic to the circle and carrying only one nondegenerate saddle 
point), the problem was settled by R. Roussarie [Rou89]. He estimated how many 
derivatives of the Abelian integral I{t) should have zero limits as t ^ t^ in order 
for the perturbation (3.2) to produce n or more limit cycles close to the separatrix 
loop on the critical level curve {H = t^}. This result yields an upper bound for 
the cyclicity of the separatrix loop in terms of the "multiplicity of the root" of the 
Abelian integral Ii at the point where Ii in fact loses its analyticity. The only 
other type of critical level curve that can occur for a generic polynomial H, if the 
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eight-shaped curve (also carrying a nondegenerate saddle). The answer in this case 
is not yet known. 

The second, apparently much more difficult problem, appears in connection 
with the natural question, How many consecutive variations Ik should be computed 
in order to guarantee that the perturbation is non-conservative? 

The inductive process described in Lemma 2, is algebraic in the following sense. 
Starting from the perturbation (3.2) given by the form u) = cu\, we construct a se- 
quence of the polynomial 1-forms lu2,(-U3, . . . which express the higher variations 
Ik{t) (3.4), assuming that all previous variations vanish identically, Ii = ■■■ = 
Ik-i = 0. Coefficients of the forms are given by algebraic (polynomial) expres- 
sions involving coefficients of the initial form u. Vanishing of their integrals implies 
an infinite number of polynomial identities between coefficients of the initial form u. 
By the Hilbert basis theorem, all these infinitely many identities are corollaries to 
only finitely many of them. Thus only finitely many steps of the inductive process 
must be performed (their number < +oo depends on the degree of the initial 
form w and, naturally, on the Hamiltonian H). If all integrals § uJk-, k = 1,. . . ,N 
arc identically zeros, then all higher variations are necessarily zeros and the entire 
family (3.2) is integrable. 

However, the problem of finding an upper bound for the number N in terms of 
iv and H is overtly open: even for the most simple case H = ^(.T^+y^) considered in 
Example 3.1 above, the answer is unknown, moreover, it constitutes the challenging 
problem on distinguishing between center and focus, posed by Poincare a century 
ago and still open. A similar though apparently more simple problem was recently 
studied in [BFY99]. 

Remark. A polynomial vector field can be integrable but not Hamiltonian: 
it is sufficient that the corresponding form fl possess an integrating factor. A 

typical example is that of Darboux integrable vector fields. Let Ai, . . . , A„ G M be 
a collection of real numbers and Hi,.. . £ IR[a;,y] polynomials. The rational 
1-form 

^^' = EA.f^ = f , F = Ht...H^'^, 

determines a conservative (integrable) singular foliation, since its null spaces {51' = 
0} are tangent to the real level curves {F{x, y) = t} C which are all closed unless 
beginning and ending at infinity. The form ft' is not polynomial (only rational), 
but the form Q = RV,' already will be, if R = gcd(i?i, . . . ,Hn) € lR[a;, t/] is the 
common divisor of all Hi . 

Appearance of limit cycles in polynomial perturbations of the form 

fl + eu = 0, ui=pdx + qdy, p, g e M[a;, y], 

is determined by the same mechanisms as in the Hamiltonian case, in particular, 
the first variation of the Poincare displacement map measured in units of In = 
Xi In Hi, is equal to the integral 

Ht)=<f ^T^W 

JF=t R{x, y) 

of the rational 1-form R~^lo along closed level curves of the transcendental function 
F. 
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As was noted on several occasions by V. Arnold, the problem on Abelian in- 
tegrals should be posed also for such "pscudo- Abelian" integrals. Nevertheless, 
there are practically no results of general nature pertinent to this problem, even 
in the simplest cases n = 2 and n = 3. In particular, all results described below 
(derivation of the Picard-Fuchs system, monodromy properties etc) fail for this 
generalized class of integrals. 

4. Finiteness problems for Abelian integrals: tangential Hilbert 

problem 

Despite its ambivalent nature, the connection between limit cycles of polyno- 
mial vector fields and isolated zeros of Abelian integrals justifies formulation of 
several finiteness problems for the latter. As suggested above, we refer to this prob- 
lem as the tangential Hilbert problem, distinguishing between several finiteness 
types according to the paradigm laid out before. 

4.1. Individual finiteness. For any fixed combination of the Hamiltonian H 
and the 1-form co the Abelian integral appearing in (3.3) as a function of the real 
variable t can be shown to be real analytic except for finitely many values of t. Ac- 
cording to §1, it may have only finitely many isolated zeros unless they accumulate 
to one of these exceptional values. The individual finiteness problem for Abelian in- 
tegrals is to prove that such accumulation is in fact impossible. This assertion turns 
out to be a simple corollary to a general theorem describing ramification of Abelian 
integrals after analytic continuation into the complex domain. The key step in the 
proof of individual finiteness for Abelian integrals is the following representation. 

Lemma 3. Any Abelian integral near any exceptional point admits a con- 
verging representation of the form 

I{t)=J2hkr{t){t-t.yin''-\t-t^) (4.1) 

r.k 

with finitely many terms, where all exponents r are rational numbers, the powers k 
take only finitely many natural values, and all functions hkr are real analytic at the 
point t = tf 

This Lemma is explained in §5. As soon as the representation (4.1) is estab- 
lished, the rest of the proof is as described in Example 1.2. It remains only to notice 
that the convergence of the expansion (4.1) is crucial: if all terms in this expansion 
are zeros, then the integral itself is identically zero and hence has no isolated roots 
at all (compare with the Remark in §2.3). 

4.2. Existential finiteness. Abelian integrals depend on the choice of the 
Hamiltonian and the form; in order to make the parameter space finite-dimensional, 
one has to restrict their degrees. Then the coefficients of the Hamiltonian H and 
the form lo become the natural parameters of the problem. 

The parameter space is intrinsically compact: indeed, replacing H and uo by 
CiH and C2W with Ci,C2 nonzero constants, clearly does not affect the number of 
isolated roots of the corresponding integral. Hence the parameter spaces are in fact 
projective spaces of appropriate dimensions. 

Fix any two integer numbers n and d and consider all Hamiltonians of degree 
< n and all polynomial forms of degree ^ d. The existence of uniform upper 
bounds for the number of isolated zeros of the Abelian integrals subject to the 
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above restrictions on the degrees, was proved by Varchenko [Var84] and Khovanskii 
[Kho84]. 

Theorem 3 (A. Varchenko- A. Khovanskii, 1984). For any n,d < +oo, the 
number of ovals 7 C {H = const} yielding isolated zeros to the Abelian integral 
§^ijj, is hounded by a constant N = N{n,d) uniformly over all Hamiltonians of 
degree < n and all polynomial 1-forms of degree < d. 

One key tool in their proof is again the same Lemma 3 properly generalized to 
cover the parametric construction. The second principal ingredient is the Pfaffian 
elimination technique [Kho91]. The latter is a method allowing to reduce the ques- 
tion on zeros of functions involving real powers, logarithms, arctangents and other 
functions that can be defined using Pfaffian equations with algebraic right hand 
sides, to the question on zeros of certain auxiliary systems of equations involving 
only algebraic functions. Application of the generalized Lemma 3 allows to derive 
the existential finiteness theorem for Abelian integrals, using in the standard way 
the compactness and localization arguments, from the similar existential finiteness 
assertion for converging multivariate expressions of the form (4.1). The latter as- 
sertion can in turn be derived, using the PfafEan elimination technique, from the 
existential finiteness theorem for analytic families of functions as introduced in §1. 
This last case was already settled by Theorem 1, and this completes the proof. 

However, since both Theorem 1 and its corollary. Theorem 2 do not give any 
information on the number of isolated roots of functions. Theorem 3 inherits this 
quality as a purely existential statement. 

4.3. Constructive finiteness. The simplest (trivial) case when constructive 
bounds for the tangential Hilbert problem can be easily produced, is the case of qua- 
dratic Hamiltonians. In this case all Abelian integrals can be explicitly computed, 
being actually polynomial functions of t. 

The constructive finite cyclicity problem for Abelian integrals is also a relatively 
simple assertion: the maximal multiplicity of an isolated zero of an Abelian integral 
admits an upper bound in terms of degw and degH found by P. Mardesic in 
[Mar91]. 

Historically the first nontrivial problem concerning zeros of Abelian integrals 

appeared in connection with what later became known as the Takens Bogdanov 
bifurcation. In order to prove that no more than one limit cycle appears by defor- 
mation of the generic cuspidal singular point on the plane, one has to verify that for 
any 1-form ui = [a + (3x)y dx its integral over the closed ovals of the Hamiltonian 
H{x, y) = y^ -\-x^ — Zx has no more than one isolated zero. R. Bogdanov in [Bog76] 
proved this claim, achieving one of the first results on Abelian integrals. Later his 
proof was considerably simplified by Ilyashenko in the paper [Ily78] which intro- 
duced some very important tools, among them the idea of complexification of the 
Abelian integrals. 

There are numerous studies treating other low-degree cases that mostly ap- 
peared in connection with bifurcations of polynomial vector fields of low degrees. 
Starting from the paper by Bogdanov, this direction was pursued, among other, by 
F. Dumortier, A. GasuU, L. Gavrilov, F. Girard, E. Horozov, I. Iliev, Yu. Ilyashenko, 
A. Jebrane, B. Li, C. Li, J. Llibre, P. Mardesic, G. Petrov, R. Roussarie, C. 
Rousseau, Z. Zhang, Y. Zhao, H. Zoladek, to mention only some names and the 



16 



SERGEI YAKOVENKO 



Hamiltonian H 


Degree of the form 


Low degree 

(degw < 2) 


Arbitrary 
d = deg uj 
(asymptotic) 


Arbitrary 
degree 
(constructive) 


Quadratic: 

H = x^+y^ 


Integrals are polynomial functions of t 


Elliptic 

H = y'^+x^ -3x 


< 1 root 
Petrov [Pet86] 




nonoscillation 

theorem 
Petrov [Pet88] 


General 
cubic H 


numerous results 

in particular, 
Gavrilov [GavOl] 




5(d+2) 
Horozov-Iliev 
[HI98] 


Hyperelliptic 

H = y^ + P^{x) 




Petrov [Pet90] 


tower function, 
Novikov— Yako- 
venko [NY99a] 


General 
nth degree 
polynomial 




C{n)d + 0{1), 
Khovanskii 
[Kho95] 


??? 



Table 2. Constructive tangential Hilbert 16th problem: partial synopsis 



Dashes in the cells mark irrelevant (artificial) problems. Blank cells indicate that 

the corresponding problem possesses no specific (more accurate) solution other 
than the one implied by stronger versions appearing to the right or down in the 
Table. 

most recent works. Among these accurate bounds, the following result is remarkable 
by its succinct formulation and difficult proof. 

Theorem 4 (L. Gavrilov-E.Horozov-I.Iliev [GavOl]). For any cubic Hamilt- 
onian H with four distinct critical values, and any quadratic 1-form w, the corre- 
sponding Abelian integral has no more than two isolated roots. 

Moreover, in the perturbation (3.2) no more than 2 limit cycles may appear, 
including those born from separatrix polygons. 

This settles the localized version of the Hilbert 16th problem for quadratic 
vector fields arbitrary close to Hamiltonian quadratic vector fields with the specified 
Hamiltonians (cubic with 4 distinct critical values). The bound is accurate. 

Another remarkable result is due to G. Petrov [Pet88] who studied completely 
the elliptic case H{x,y) = y^ + Pzix), where P3 G R[x] is a cubic univariate 
polynomial: by affine transformations one can always reduce P3 to the form P^ix) = 
x^ — 3x. Petrov proved that the Abelian integrals of forms of arbitrary degree form a 
non-oscillating (in other languages disconjugate, or Chebyshev) family: the number 
of isolated roots never becomes equal or exceeds the dimension of the linear space 
of all such integrals. 

Remark. One can easily verify that any linear space spanned by d linear inde- 
pendent analytic functions fi(t), . . . , fd{t) contains a nontrivial linear combination 
Ci/i H h Cdfd exhibiting a root of multiplicity — 1 at any preassigned point. 
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Methods introduced by Petrov were further elaborated and refined. The follow- 
ing bound obtained by E. Horozov and I. Ilicv, though not sharp, covers the case 
of Abelian integrals of arbitrary polynomial forms over level curves of any cubic 
Hamiltonians. 

Theorem 5 ([HI98]). For any cubic Hamiltonian H{x, y) the Abelian integral 
of a form of degree d cannot have more than 5{d+2) isolated zeros. 

By similar methods the quartic Hamiltonians (of degree 4) with elliptic level 
curves were studied in [GJ98]. The case of quartic elliptic Hamiltonians perturbed 
within the Lienard equation was a subject of four recent preprints of C. Li and 
F. Dumortier. 

4.4. Asymptotic bounds. Results of the different type, asymptotic bounds, 
take advantage of the asymmetry of the roles played by 1-forms and Hamiltonians. 
Unlike the original Hilbert 16th problem, in which the coefficients of the polynomials 
P, Q all enjoy equal rights as parameters, the roles of the Hamiltonian H and the 
form Lo in the definition of the Abelian integrals are fairly different (e.g., the integral 
depends on co linearly, whereas even small variations of H can result in drastic 
changes of the domain of definition). Thus it makes complete sense to separate 
these two parameters and study first the dependence on w, treating H as fixed 
("individual"). 

The corresponding "semiconstructive" problem was addressed in a number of 
recent publications. First, very excessive (double exponential in d) upper bounds 
were obtained in [IY95a] and almost immediately improved to simple exponential 
expression in [NY95] by Ilyashenko, Novikov and the author. The ultimate result 
in this direction, a linear upper bound of the form Ci{n)-d+C2{n), was obtained by 
Khovanskii [Kho95] using some ideas developed earlier by G. Petrov [Pet90]. Here 
the constant Ci(n) is absolutely explicit (e.g., does not exceed a double exponential 
of n), while the second constant C2(n) is purely existential though uniform over all 
Hamiltonians of degree < n. 

4.5. Hyperelliptic integrals. The only particular case covering Hamilton- 
ians of arbitrarily high degrees, for which constructive solution of the tangential 
Hilbert problem is known, is that of hyperelliptic Hamiltonians, 

H{x,y) =y'^ +Pn+i{x), P„+i = a;"+^ + aia;""^ + • • • + a„_ia; + o„. (4.2) 

Singular points of the corresponding Hamiltonian system correspond to critical 
points of the polynomial P, called potential. The reason why the class of hyperel- 
liptic polynomials is especially simple, one can vaguely attribute to the fact that 
from many points of view, hyperelliptic polynomials behave like their univariate 
potentials. In particular, this concerns topology of the bundles defined by complex- 
ification of H. 

Under the addition assumption (believed to be technical, though it occurs in- 
dependently and persistently in several related problems) that all critical points of 
the potential P are real, the tangential Hilbert problem turns out to be construc- 
tively solvable. More precisely, as shown by D. Novikov and the author [NY99a], 
there exists an algorithm defining an elementary function C{n,d) of two integer 
arguments n and d, such that for any form of degree ^ d and any hyperelliptic 
Hamiltonian n of degree < n having only real critical points, the corresponding 
hyperelliptic Abelian integral has no more than C{n, d) isolated roots. 
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The algorithm involves several nested inductive constructions, resulting in an 
extremely excessive bound: it is given by a tower function (an iterated exponential) 
of height greater than 5 but probably smaller than 10. 

The proof is based on the fact that Abelian integrals satisfy a system of first 
order linear ordinary differential equations with rational coefficients, called the 
Picard Fuchs system. 

4.6. Quantitative theory of ordinary differential equations as a tool 
for constructive solution of the tangential Hilbert problem. The method 

used for constructive solution of the tangential Hilbert problem in the hypcrcUiptic 
case, is fairly general. In [Yak99] one can find an introduction to the general theory 
allowing to investigate zeros of functions defined by ordinary differential equations 
with polynomial and rational coefficients. Basics of this theory were developed in a 
series of joint works by D. Novikov and the author; they are briefly recalled below. 

The principal goal of these lecture notes is twofold. First, we show that the 
Picard Fuchs system of differential equations can be explicitly written down for 
an arbitrary generic Hamiltonian (not necessarily a hyperelliptic one), at least at 
the price of certain redundancy. Some additional information may be extracted 
from the explicit derivation procedure. In particular, we show that for almost any 
Hamiltonian H (chosen for convenience of degree n + 1): 

1. there exist /U different Abelian integrals Ik = § cok of monomial 1-forms 
LVi, . . . ,LOfj,, n = n^, satisfying together a Fuchsian system of linear ordinary 
differential equations with rational coefficients of the form 

|.,«,, m^t^. ^^^^ 

xeC, Aj- e Mat^x/.(C), t e C,ti,...,i^ G C, 
such that 

2. integral I{t) of any other form u) can be represented as a linear combination 
of the integrals Ji, . . . , J^, 

lit) = gi {t)h {t) + --- + q^ {t)I^ it) , (4.4) 

with polynomial coefficients Qj S C[f], degg^ ^ degw/degiJ. 
The residue matrices Aj can be sufficiently completely described, in particular, 
upper bounds on their norms can be placed. 

The second goal is to introduce several results on the number of zeros of func- 
tions defined by systems of ordinary linear equations with polynomial and rational 
coefficients, and polynomial combinations of such functions. We start with the 
simplest case of one nth order linear equation and show that isolated zeros of its 
solutions can be described in terms of the magnitude of coefficients of this equation. 
In the simplest real case this is a classical theorem by de la Vallce Poussin [dlVP29], 
which we generalize for the complex analytic context. Afterwards we study systems 
of polynomial ODE's in the real and complex space of arbitrary dimension; here for 
the first time appear explicit bounds given by tower functions. Finally, we discuss 
the case of Fuchsian systems of the form (4.3) and show that under certain natural 
restrictions on the monodromy group of this equation, a global upper bound on the 
number of zeros can be given in terms of norms of the residue matrices Aj . 
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One can hope that combination of these two techniques ultimately would allow 

to construct an explicit bound for the tangential Hilbert problem in the general 
case, filling the right bottom corner of the Table 2. 

Lecture 2. Abelian integrals and differential equations 

In this lecture wc explain the connection between Abelian integrals and linear 
ordinary differential equations. 

5. Complexification of Abelian integrals: topological approach 

5.1. General scheme. We recall here the basic construction of complexifica- 
tion of Abelian integrals. All details can be found in various textbooks, [AGV88] 
being the principal source. 

A bivariate polynomial H £ C[x, y] defines a map from to with preimages 
of points being affine algebraic curves. It turns out that the map H is a, topological 
bundle over a complement to finitely many points E = {ti, . . . ,tr} C C. This 
allows to identify in a canonical way the homology groups of all fibers H~^{t) 
sufficiently close to H~^{t^), which in turn allows to extend integrals over 1-cycles 
of polynomial 1-forins as complex analytic functions of t in U, ramified over the 
singular locus S. Geometrically this can be expressed as introducing a locally flat 
connexion on the (co)homological bundles over the punctured sphere. 

Prom the same topological arguments it follows that for any polynomial 1- 
form CO the linear space generated by integrals ui over any family Si{t), . . . , S^^{t) 
of cycles forming a basis in the first homology group of the fiber Xt — H^^{t), 
is invariant by analjd;ic continuation along the loops avoiding the locus E. The 
monodromy group consisting of all automorphisms of this space occurring as the 
result of continuation over all loops, is independent of the choice of the form uj and 
depends only on the Hamiltonian H. 

These fairly general topological considerations already imply that the Abelian 
integrals form a finitely generated module over the ring of polynomial functions 
of t, and generators of this module satisfy a system of first order linear ordinary 
differential equations with rational (in t) coefficients, called Picard-Fuchs equations. 

However, in order to apply methods described in subsequent sections, it is 
necessary to obtain quantitative characteristics of these equations, in particular 
their dimension, degree and the magnitude of coefficients. Part of this information 
(e.g., the dimension) can be easily achieved from the above construction. To obtain 
upper bounds on the degrees, some other more refined considerations are required, 
but it is practically impossible to derive bounds for the coeflicients using only 
topological arguments as above. 

The current section contains a brief exposition of well-known facts leading to 
derivation of Picard-Fuchs equations and representation of the space of Abelian 
integrals as the Picard-Vessiot extension. In the next section we suggest an alter- 
native approach based on elementary algebraic consideration, that allows to derive 
explicitly the Picard-Fuchs system at the price of certain redundancy. 

5.2. Topological bundles defined by proper maps. Let /: M ^ A'^ be a 

smooth map between two manifolds. Recall that a point & e A'' is a regular value 
for /, if the rank of the differential /* : T^M ^ T^N is maximal (equal to dim A/') at 
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all points a; e X(, of the preimage. Complement to the set of regular values consists 
of critical values and is denoted crit /. 

Lemma 4. If b is a regular value of a proper map f : M ^ N , then there exists 
a neighborhood U 5 b such that all preimages Xy = f~^{y) C M are diffeomorphic 
to Xb = f-'^{b) so that f-^{U) c^XbxU. 

Proof. Consider an arbitrary vector vq G T^N and embed it into a vector field 

V on N . Wc claim that in a sufficiently small neighborhood of the preimage X^ 
one can construct a smooth vector field w such that f^w = v, that is, w and v are 
/-related. Such a field obviously exists near each point x G Xh, since /* is surjective 
(and takes the form of a parallel projection in suitably chosen local coordinates by 
virtue of the theorem on rank). Now it remains to choose a finite covering of Xi, 
by these neighborhoods and patch together the corresponding vector fields, using 
the appropriate partition of unity. 

To conclude the proof, notice that the (local) flows of v and w are conjugate 
by / (by construction), hence the local flow of w, defined in some neighborhood 
of Xf,, takes the latter preimage into preimage of the corresponding point y on the 
flow curve of b. Since the initial vector vq can be chosen pointing to any direction, 
this proves that all sufficiently close preimages are diffeomorphic to each other. □ 

Remark. The diffeomorphism between close preimages is not canonically de- 
fined, but any two such diffeomorphisms are homotopically equivalent, since the 
vector field w from the proof of the Lemma is homotopically unique. 

Corollary 1. IfM is compact, then f is a topological bundle over N\cvitf. 

□ 

Remark. Assertion of this lemma for the case N = M. is the fundamental 
principle of the Morse theory. 

5.3. Homology bundle and flat connexion on it. A locally trivial topo- 
logical bundle f:M^N defines in a canonical way the associated homology 
bundle over N with the fibers being the homology groups Hi{Xy,Z) (we will be 
only interested in the case i = 1) together with a flat connexion on this bimdle. 
This means that any 1-cycle on any particular fiber Xa can be transported to any 
other fiber X;, along any path 7 connecting a and b in N. Flatness means that 
continuation along any sufficiently small loop beginning and ending at a, returns 
any 1-cycle to its initial position. However, transport along "long" loops (not con- 
tractible in N) can result in a nontrivial transformation of the homology, called the 
monodromy transforrn,ation. 

Fix a point a G N and choose a basis of 1-cycles in the group 

H\{Xa,7j), where fi is the rank of the homology group. Then the monodromy 
transformation corresponding to any loop 7 from the fundamental group tti {N^ a) 
can be described by the corresponding monodromy matrix M = e ^l&t^xni'^) 
with integral entries m,j : 

A-y^j = Jjmy , i = l,...,/x. (5-1) 

In the matrix form the monodromy transformation acts on the row vector 5 = 
((5i, . . . , (5ju) as multiplication by the matrix from the right, A^(5 = 8 ■ My. 
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The integer-valued matrix is obviously invertible. Moreover, the corre- 
spondence 7 1-^ My is an (anti) representation of the fundamental group Tri{N,a) 
in GL(n, Z). This implies, in particular, that detM^ = ±1 for all loops 7. 

The flat connexion on the homology bundle defines a connexion on the co- 
homology bundle which makes it possible to compute the covariant derivative of 
1-forms: any smooth 1-form a; on M can be restricted on any fiber Xa and in- 
tegrated along a continuous (horizontal, locally constant) family of 1-cycles 6{y), 
resulting in a function I{y). The result of continuation of this function along closed 
paths in N is completely determined by the monodromy group and does not de- 
pend on the choice of the form u). More precisely, integrating the form u over each 
of the 1-cycles S{t) = {Si{y), . . . , (5^(y)), we obtain a tuple of continuous functions 
Ij{y) = §^^y^LU, j = 1, . . . , /Lt, which after continuation along a path 7 G tti{N, a) 
undergo the transformation 

A^I = I.M^, J = (/i, . . . , I^), Ij = Ij{y) (5.2) 

with the same matrices independently of the form u. This basic fact lies in the 

core of the topological theory outlined below. 

5.4. Topological bundles defined by polynomial maps. Our goal is to 
apply the previous construction to the polynomial map _ff : — ^ C considered as 
a smooth map between smooth manifolds. Since compactness of the preimages is 
crucial for these arguments, we need to compactify the domain (and the range) of 
H. 

In contrast with the one-dimensional case, it is in general impossible to extend 
iJ as a map between the natural compactifications CP^ and CP^ respectively, 
since on the infinite line CP^ C CP^ one has several points of indeterminacy: they 
occur at the intersections between compactified preimages H~^{t) and the infinite 
line. Algebraically this can be seen as the indeterminacy of the rational expression 
H{l/z, y/z) = P{y, z)/z'^, d = degH, at the points where P{y, 0) = (at all other 
points of the infinite line {z = 0} one can assign the value -ff = 00 to this ratio). 

The problem can be resolved by blowing up these indeterminacy points, in the 
same way as blowing up the origin allows to assign values from CP^ = C U {00} 
to the rational expression R{x,y) = y/x that is initially indeterminate. After a 
series of blow-ups at indeterminacy points one arrives at a compact two-dimensional 
complex manifold M and a map (still denoted by H) from M to CP^, called 
determination of the initial polynomial map. 

Now one can apply Lemma 4. Since the determination H: M ^ CP^ is alge- 
braic, it has only a finite number of critical values Sh = {ti, - ■ ■ ,ts} C CP^ and we 
conclude that it defines a topological bundle over the complement of these excep- 
tional values. Any 1-cycle S{a) C H~^{a) on a fiber of this bundle can be embedded 
into a continuous horizontal family 5{t) of 1-cycles, ramified over the exceptional 
locus. A polynomial 1-form lo extends as a meromorphic 1-form on CP^ with the 
polar divisor CP^^ that after the blowing up becomes an algebraic hypersurface 
D in M. Let E be the union of Sh and the critical values of the projection H 
restricted on D (including the images of the non-smooth points of D). 

The puUback of a; on M can be integrated along the family 5{t): by the Cauchy- 
Stokes theorem, the result depends only on the homology class of the cycle. It can 
be easily shown that the result of this integration is an analytic function of t, 
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ramified over E. Its monodromy (the result of analytic continuation along closed 
loops avoiding the exceptional locus S) is as before determined only by H. 
This construction proves the following result. 

Proposition 1 (cf. with [AGV88]). For any polynomial if: C, any 

1- cycle 5 on a nonsingular level curve = {H = t^,} C and any polynomial 

1-formuj the Abelian integral I {t) = ^^uj can he extended as an analytic multivalued 
function ramified over a finite number of points depending only on H . □ 

Behavior of the integrals near the ramification locus is relatively tame: it can 
be seen that any integral can grow no faster than polynomially in \t — tj\~^ as t 
tends to some tj G S remaining in any sector with the vertex at tj. 

5.5. Picard— Fuchs system. Let as before 5{t) = {Si{t), . . . , (^^(f)) be a con- 
tinuous family of 1-cycles forming a basis (frame) of the first homology group of 
the respective fibers H~^{t), arranged as a row vector. One can show that there 
exist /i polynomial 1-forms wi , . . . , such that the period matrix 



(5.3) 



is not identically degenerate, AciX{t) ^ 0. From (5.2) it follows that 

A,^X(t) = X{t)M^, V7 e 7ri(Cpi \ S,a). (5.4) 

Differentiating the identity (5.4), we see that the derivative X{t) has the same 
monodromy (i.e., X is multiplied by the same matrix factors M^). Therefore the 
"logarithmic derivative" A{t) = X(t) ■ X^'^(t) is single-valued (invariant by all 
monodromy transformations) meromorphic matrix function having poles of finite 
order at the points of S and eventually at the points of degeneracy of X{-): 

A^A{t) = X{t)M^ ■ M-^X-'^{t) = X{t) ■ X-'^{t) = A{t) 

for any loop 7 e 7ri(C \ i7, t*). Prom this we conclude that A{t) is a rational 

matrix function while the period matrix X{t) is a fundamental matrix solution 
to the system of linear ordinary differential equations with rational coefficients on 
CP\ 

X = A{t)X, A{-) e Mat^x^(C(t)). (5.5) 

The common name for various such systems satisfied by Abelian integrals, is the 
Picard-Fuchs system (or Picard-Fuchs equation). 

Integrals of any other form oj can be expressed as linear combinations of in- 
tegrals of the framing forms cj,;, i = 1, . . . with coefficients from the field C{t,) 
of rational functions. Indeed, multiplying the row vector I{t) = {Ii{t), . . . ,/^(i)), 
Ij — §^ IV, by X~^{t), we obtain (for the same reasons as above) a single-valued 
hence rational row- vector function r{t) = {ri{t), . . . ,r^(t)), that is. 



Js(t) ^ J5(t) 



r, e C(i). (5.6) 



for any cycle 5{t) continuously depending on t. Note that the space of functions 
representable as (5.6), is closed by derivation. 
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We summarize this as follows. Recall that a Picard-Vessiot extension is a 
differential field of analytic multivalued functions obtained by adjoining to the field 
C{t) all components of a fundamental matrix solution of a system of linear ordinary 
differential equations with rational coefficients. 

Theorem 6. Abelian integrals belong to a Picard-Vessiot extension for some 
system of linear ordinary differential equations with rational coefficients. 

Later we discuss this formulation and relevant issues in more details. However, 
it is important to stress here that neither entries of the rational matrix A(t) nor the 
rational coefficients rj{t) can be computed explicitly from the above construction 
without additional considerations. 

6. Hamiltonians transversal to infinity 

6.1. Definition. For an arbitrary Hamiltonian H, even location of the rami- 
fication points tj is difficult to describe without effectively resolving all the indeter- 
minacy points at infinity. However, under some natural and generic assumptions 
one may guarantee that no new critical points will appear after compactification 
and blowing up the indeterminacy points on the infinite line CP^ C CP^. This 
will immediately imply that the set S of exceptional values must be a subset of 
crit H, the set of critical values corresponding to finite critical points in only. 

Definition 5. A polynomial H e C[x,y] is said to be transversal to infinity, 
if its principal homogeneous part L = J2i+j=d^ij^^y'' ' ~ degi?, factors as a 
product of pairwise different linear forms. 

Equivalent conditions follow. 

f . The principal homogeneous part has an isolated critical point at the origin; 

2. The partial derivatives ^ are mutually prime; 

3. H has exactly = {degH — 1)^ critical points in if counted with multi- 
plicities; 

4. Each level curve {H = t} intersects transversely the infinite line 'CP^ C CP^ 
after projective compactification. 

6.2. Topology of polynomials transversal to infinity. We prove now that 
for H transversal to infinity, E = crit H. 

Proposition 2. A polynomial H: ^ transversal to infinity is a topo- 
logical bundle over the set crit of critical values of H . 

Instead of proving this by resolving the indeterminacy points at infinity, one 
may modiiy the proof of Lemma 4 sketched above, and construct the vector field 
w near an arbitrary infinite point p G CP^ on CP^ with the property H^w = ^. 

Proof, ft will be shown below in §7.1 that if H of degree n + 1 is transversal 
to infinity, then there exist two polynomials a,b G C[x,y] of degree n — 1 such that 

a— + b— =x^" ' + •••, 
ox ay 

where the dots stand for a bivariate polynomial of degree ^ 2n — 2. One can easily 
check that the rational vector field 

^ ^ a{x,y) d ^ b{x,y) d 
2;2n-i -)-... Qx a;^"~i + ■ ■ ■ dy' 
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i?-related to the field in tlie cliart (l/x, y/x) is regular (smooth) on the infinite 
line l/x = Q (more precisely, on its afRne part covered by this chart). The other 
affine part is covered by the field that is obtained in a similar way from solution of 
the equation a|f = 2/^""^ + • • • • 

Thus near each point of the compactified level curve Xa C CP^ one has a 
smooth vector field if-related to ^ in the finite part C CP^ (in particular, this 
implies that this field vanishes at all points of indeterminacy of H on CPj^). The 
rest of the proof is the same as in Lemma 4. □ 

6.3. Module of Abelian integrals, Gavrilov and Novikov theorems. 

For Hamiltonians transversal to infinity, the constructions of §5 can be further 
refined. In particular, the choice of the forms wi, . . . can be made much more 
explicit. 

Lemma 5 (L. Gavrilov [Gav98]). Let H he a Hamiltonian of degree n + 1 
transversal to infinity, with distinct critical values /x = n^. Then one 

can choose v? monomial 1-forms Wi, . . . ,a;^ of degrees ^ 2n so that the respective 
period matrix has the determinant detX{t) = c{t — ti) ■ ■ ■ {t — t^) with c^O. 

In this assertion we use the convention on degrees of polynomial A -forms for- 
mulated in §7.1. As a corollary, one can derive the following result that refines the 
assertion of Theorem 6. 

Corollary 2. Abelian integral of a 1-form to of degree d can he represented 

as 

f w = V'pj(i) f ujj, 

J5(t) ^ Js(t) (6.1) 

Pj{t) e C[i], degWj + AegH ■ degpj < degw. 

In other words, Abelian integrals constitute a module over the ring C[f] that is 
generated by integrals of the basic forms ujj. 

The constant c from Lemma 5 depends on the choice of the monomial forms 
and the Hamiltonian. Its value was explicitly computed by A. Glutsuk [GluOO] 
following some ideas of Yu. llyashenko, and a simple elementary proof of the in- 
equality c 7^ for an appropriate choice of the monomial forms was obtained by 
D. Novikov [NovOla]. In the same paper [NovOla] it is proved, using some of 
the methods described below, that the period matrix satisfies a system of linear 
ordinary differential equations 

with a matrix polynomial P(t) of degree /i, in general having Fuchsian singularities 
at all points tj of the ramification locus S, but a non-Fuchsian singularity &tt = oo. 

As yet another corollary, one can derive a Picard-Puchs system for the period 
matrix. Let LOj be as in Corollary 2. Consider the closed 2-forms dH A uij and let 
flj be any polynomial primitives satisfying the conditions 

dflj = dH A LOj, j — 1, . . . , 1^. 
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Each 0, can be expanded as in Corollary 2, yielding an identity between the inte- 
grals, 

i=i 

valid for any choice of a continuous family of cycles of integration 5{t). On the 
other hand, as will be shown in §7.2, the derivative of each integral § is exactly 
the integral § tOi. Differentiating the above identities, we arrive to the matrix 
differential equation 

X = PX + PX, P = P{t) = \\p^,{t)\\ G Mat^xM(C[t]) 

with a matrix polynomial P{t) of some known degree. If required, this identity can 
be resolved to the form (5.5). 

6.4. Commentaries. The proof of Lemma 5 and Corollary 2 is based on a 

finer than before analysis of topology of the bundle H for polynomials transversal 
to infinity. In particular, in assumptions of the Lemma, one can choose a special 
framing of the homology bundle by vanishing cycles Sj{t), represented by loops 
on the prcimagc Xt that shrinks to a point when t tj (a special precaution is 
required to avoid problems with multivaluedness) . For such choice of the cycles, 
the period matrix X{t) must have a vanishing column at each of the points tj, for 
any collection of the framing forms. Next, in this case the determinant dct X is a 
single-valued function that therefore must be a polynomial (being locally bounded 
everywhere on C). Its growth as t ^ oo depends on the degrees of the framing 
forms, since the "size" of the cycles Sj{t) grows in a known way (depending only 
on H). 

All this implies that the determinant of any period matrix detX{t) is a poly- 
nomial in t divisible by ]~[j=i(* ^ ^j)- The coefficients Pj(t) of representation (6.1) 
can be found by solving a system of linear algebraic equations. When solved by the 
Cramer rule, this system yields pj (t) as a ratio of two determinants, the determi- 
nant of the period matrix for the basic forms tJi, . . . in the denominator, and 
that of a similar period matrix for the collection of 1-forms with ujj replaced by u) 
in the numerator. By Lemma 5, these ratios are polynomials in t and their degrees 
can be easily majorized in terms oi d = degw. 

6.5. Multiplicity of the roots of Abelian integrals: constructive finite 
cyclicity. A similar construction allows to majorize the order of a zero of any 
Abelian integral, at least under the assumption that H is transversal to infinity 

and has only Morse critical points [MarOl]. Instead of the period matrix X, 
consider the matrix function J{t) whose entry Jij{t) is the {i — l)-st derivative of 
the integral Ij{t) = §s -(t) 

For the same reasons as before, its determinant w{t) = det J(i), the Wronskian 
of the integrals Ii,. . . ,1^, is a single- valued hence rational function of t. Its poles 
may occur only at the points tj and t = oo. The assumption on finite singular 
points implies that all integrals Ij{t) have at worst logarithmic growth near each 
tj, and this growth rate allows for differentiation so that /j* ^\t) grows no faster 
than \t — tj\^~^ ast tj without spiraling. These estimates imply an upper bound 
on the total order of all poles of w at all finite points. The growth rate of w as 
f — > oo depends on deg co and can be easily estimated. This gives an upper bound 
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on the degree v = degw(f) of the rational function (the total number of its poles on 
CP^ including those at infinity). This degree is obviously an upper bound for the 
order of any nontrivial zero of w. Since the order of the Wronskian w is equal to 
/U, the number v + ^—1 is an upper bound for the order of any root of any integral 
Ij{t) at any point t^tj. 

6.6. Reservations. Despite more detailed constructions and more accurate 
considerations, the approach based only on topological ideas cannot provide many 
important data. For example, even the "constant" c from Lemma 5 depends in 
a rather nontrivial way on both H and the choice of the framing forms ujj, see 
[GluOO] and [NovOla]. Among other things, this means practical impossibility of 
majorizing the polynomial coefficients pj in (6.1). The same refers to the derivation 
of the Picard Fuchs system: after an accurate computation, it can be reduced to 
the form determined by three constant matrices Pi, P2, P3, 

{Po + tPi)X = P2X, Pi G Mat^x^(C) 

but no bounds on the norms of these matrices (or their inverses) can be obtained 
except for some especially simple cases. 

In short, it is combinatorial parameters like degrees, dimensions, ranks and so 
on that can be more or less easily derived from even the most explicit topological 
constructions. In contrast, all magnitude-like parameters (norms, absolute values, 
diameters of point sets etc.), require additional arguments. Some of them can be 
obtained using rather elementary algebraic considerations. 

7. Elementary derivation of the Picard Fuchs system 

In this section we derive yet another Picard-Fuchs system by explicit linear al- 
gebraic considerations. The advantage of this approach (besides its transparency) 
is that it allows to bound explicitly the magnitude of coefficients of the system. In 
addition, the Picard-Fuchs system obtained this way possesses a nice hypergeomet- 
ric form, exhibiting only Fuchsian singularities (though this fact was not yet fully 
exploited). However, this transparency and explicitness is achieved at the price of 
a certain redundancy: the dimension of the system obtained this way, is two times 
bigger than the minimal possible one. 

7.1. Division by the gradient ideal. Division with remainder by an ideal 

can be expressed in the language of polynomial differential forms. In what follows we 
consider A;-forms A'^, k = 0,1, 2, with polynomial coefficients on the plane C^. The 
degree of a fc-form by definition is the maximum of degrees of its coefficients, plus 
k. Under such convention, deg(^Ar7) ^ deg^-t-degry for all admissible combination 
of ranks of ^ and 77 between and 2, and also 

degdw ^ degoj Vw e A*', 

for any rank k = 0, 1. The linear space of forms of rank k and degree d will be 

denoted by A^, and we denote A|^ = + A'[ -\ h A^. 

Consider co = adx + bdy £ with only isolated singularities. This implies 
that the ideal (a, 6) generated by the polynomials a,b G C[x,y] ~ A° has a finite 
codimension fj, in C[x,y], that is, there exist fj, polynomials ri,. . . ,rfj, G C[x,y] such 
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that any polynomial q from this ring can be represented as 

q = av - bu + ^^Ciri, u,v € C[x,y], Ci,...,c^gC. 
1 

Introducing 2-forms Q ~ qdx A dy, Ri = r-idx A dy, i = 1, . . . , /i. and 1-form 
r] = udx + vdy the above identity can be rewritten as division with remainder, 

VfieA^ 3r]€A^: n = wAr] + R, R = ^CiRi&K^. (7.1) 

i=l 

In particular, any 2-form O can be divided with remainder by the differential 
dH of any polynomial H ^ hP . Note that the space of remainders may be arbitrarily 
enlarged if necessary: the uniqueness of the division (7.1) will be lost then, but in 
exchange one may get better norms of the ratio and remainder. 

In general the procedure of division can be very delicate. However, if H is 
transversal to infinity, then one can easily describe the outcome, explicitly majoriz- 
ing the degrees of the remainder R G A'^ and the incomplete ratio t] G A^. 

Lemma 6 (sec [NYOlb]). If the polynomial H G A" of degree n+1 is transver- 
sal to infinity, then any 2-form CI can be divided by dH with the incomplete ratio rj 
of degree ^ dcgO — dcgH and the remainder R of degree ^ 2n. 

Proof. If L g A° is a homogeneous polynomial of degree n+1 without multiple 
linear factors, then the map between subspaces of homogeneous forms, 

;i = aL: A^^AL+1, r]^dL AT], (7.2) 

is an isomorphism. Indeed, in the bases consisting of all monomial forms of the 
given degrees, the matrix of -3 is the Sylvester m,a,trix whose determinant is the 
resultant of the two partial derivatives and Ly. The assumption on H implies 
that this resultant is nonzero. 

Therefore any homogeneous form of degree exactly 2n + 1 is divisible by dL. 
Any monomial form of degree greater than 2n + 1 can be represented as a monomial 
1-form of degree 2n + 1 times a monomial function and hence is divisible by dL 
without remainder with the same relation between the degrees, 

dcg n ^ 2n + 1 =^ CI = dL Arj, dcg rj ^ dcg fl — deg L. 

Applying this observation to all homogeneous components of a 2-form = fio -|- 

h ^2n + ^2n+i H \-fld = R + ^2n+i H h ^d, wc provc the assertion of the 

Lemma for homogeneous polynomials. 

To divide a 2-form f2 of degree ^ 2n -|- 1 by a nonhomogeneous differential 
dH = dL + ^, where L is the principal homogeneous part of H, deg^ < degL, we 
divide it by dL first, and then transform the result as follows, 

n = dLA'n + R^{dH-(,)AT] + R = dHAT] + n', 

fi' = i? - ^ A ?7, deg n' max(deg ft - dcg ^ + dcg L, deg R) < deg O, 

reducing division of by dH to division of another form f2' of strictly inferior degree. 
Iterating this step, we prove the Lemma in the general case. Notice that this is 
essentially the algorithm of division with remainder of univariate polynomials. □ 

Remark. Lemma 6 is an example of the redundant division. For a Hamiltonian 
H of degree n+1 transversal to infinity, the gradient ideal (ideal of 2-forms divisible 
by dH) has codimension n^. Indeed, the codimension is equal to the number of 
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critical points of H in C^, counted with their multipHcities. This latter number 
is exactly by virtue of Bczout theorem, since no critical points are allowed to 
"escape to infinity" by the assumption on H. 

On the other hand, the linear space of bivariate monomials of degree < 2n — 2 
(the space of 2-forms A^jn) is = 2n(2n — l)/2 w 2n^, roughly two times greater 
than fi. 

The irredundant analog of this theorem can be easily restored if necessary. 
However, the choice of monomial 2-forms generating the remainder, will depend on 
the principal homogeneous part L oi H. 

7.2. Gelfand Leray residue and derivative. 

Lemma 7. Let Lu,r] € be two polynomial 1-forms such that 

duj = dH A rj. 

Then for any continuous family S{t) of 1 -cycles on the level curves H~^{t), 

d 



U) = f T]. 
Js{t) J&{t) 

The proof in the real case (assuming only smoothness of the forms) can be 
achieved by integration of duj over the annulus between 5(t) and 5{t + At) and 
passing to limit as At — > 0. 

This formula allows to differentiate explicitly Abelian integrals of a form lo, 
expressing the result as an Abelian integral once again if du is divisible by dH. In 
fact, T} can be a rational 1-form having all zero residues after restriction on each 
level curve H~^{t). 

Example 7.1. Let H{x,y) ~ + y"^ and lo = ydx. Then §jj^^'jJ = -Trt 
(consider only the real values of t and use the Stokes formula for the circle positively 
oriented). Clearly, the form r] = ^y~^ dx satisfies the assumption of the Lemma, 
and indeed 



JH=t Jo ^ 



"^^ dcoss 



= — TT. 



sms 

This example helps to memorize the order of the wedge multiplication in the 
Gelfand-Leray formula. -4 

7.3. Derivation of the redundant Picard— Fuchs system. The linear space 
A^2n of possible remainders occurring in the division (7.1), is spanned by mono- 
mial forms x'''if dx A dy, r + s ^ 2n — 2. Denote its dimension (2n — \)n by u 
and choose any monomial primitives Wi, i = 1, . . . ,v, so that dwi span the quotient 
space A^2n/'^A^2n (modulo exact forms). Below we refer to Wj as the basic forms. 

Consider the 2-forms H dui G A<3„_|_j and divide them with remainder by dH: 

H duji = dH A r]i + Ri, i = l,...,v. (7.3) 

By the assertion on the degrees, dcgrji ^ dcguji + dcgH — degdH ^ dcgoji ^ 2n, 
therefore each of the forms can be represented (modulo an exact polynomial form) 
as a linear combination of the basic forms, 

V 
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Similarly, being all of degree < 2n, the remainders Ri £ A<2„ can be represented 
as linear combinations of the forms duf. 

V 

Ri = ttij du)j, aij e C. 

Let 5{t) be any continuous family of cycles. Then for any t ^ S the forms 

H du>i — aijdwj , Vi = 1, . . . , z/, 

i=i 

are all divisible by dH with the ratios being cohomologous to Ylj=i ^ij^j- Denote 

Xi{t) = (t Ui, i=l,...,u. 

J5(t) 

Note that integration of a form Hw over any cycle 5 C {H = t} yields t §^ oj, since 
H is constant on the cycle. Integrating both sides over the oval S{t) and using the 
Gelfand-Leray formula, wc arrive to the identities 

tXiit) - E = E bijXjit), (7.4) 

which means that the column vector {Xi, . . . , Xi,) of Abelian integrals satisfies the 
system of linear ordinary differential equations 

{tE-A)X = BX, XeC, e Mat^xL.(C), (7.5) 

with the constant matrices A = ||ay||, B = as parameters {E is the identity 

matrix). Writing 

1 

[tE - A)-i = — • P{t), P = E Pkt\ X{t) = det{tE - A), (7.6) 
^^^> fe=o 
where P{t) is a matrix polynomial of degree v —I, the adjugate matrix for tE — A, 
one sees immediately that (7.5) is a system of linear ordinary differential equations 

with rational coefficients. 

7.4. Hypergeometric systems. The form (7.5) is rather specific: for in- 
stance, all singular points of this system are Fuchsian (simple poles), including the 

point at t = oo. This is obvious if the spectrum of A is simple, but holds true in 
the general case as well, as follows from the explicit formula for [tE — A)~^ ior A 
in the Jordan normal form. 

It would be appropriate to remark here that the residues Aj of the Fuchsian 
system (4.2) are invariant by any conformal change of the independent variable 
t. In the case of hypergeometric systems (7.5) the point f = oo is distinguished: 
the residues Aj at all finite points tj G Spec A have rank 1 for a generic matrix 
A, whereas the rank of the residue A^ = — Aj at infinity is generically full. 
Thus the natural symmetry group of hypergeometric systems is not the full group 
of conformal automorphisms of CP^ , but rather the affine group of transformations 
t ^ at + b, a,b € C fixing the point t = oo. Making an affine transformation 
transforms the system (7.5) into the system {tE + A')X = B'X with the same 
matrix B' = B and A' = a-'^{A- bE). 
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In the subsequent sections it will be shown that in order to estimate the number 

of isolated zeros of solutions to Fuchsian systems, it is sufficient to know the norms 
of the residue matrices. As follows from the explicit inversion formula (7.6), norms 
of the residues can be bounded if the norm |m| is bounded from above and pairwise 
distances between the critical points tj are bounded from below. Of course, choosing 
a suitable afEne transformation as above, one can change the norm \\A\\, but at 
the same rate the distances between the singular points (eigenvalues of A) will be 
affected. In other words, the norm of the matrix ||A|| should be majorized relative 
to the spread of its eigenvalues. 

7.5. Explicitness and bounds. The above derivation does not involve any 
existential assertion: all constructions are completely transparent and allow for 

explicit bounds, say, on the norms of the matrices A, B from (7.5). To do this, 
we introduce the norms on the ring of polynomials, letting \\p\\ being the sum of 
absolute values of all its coefficients. This norm is multiplicative, \\pq\\ ^ \\p\\ \\q\\, 
and extends on polynomial fc-forms, remaining multiplicative with respect to the 
wedge product. The exterior derivative is a bounded operator on forms of bounded 
degrees. 

The entries of the matrices A, B appear as coefficients of linear expansion of 
a known 1-form in the chosen basis. But since the basic forms uji are monomial 
with coefficients equal to 1, to majorize these entries it is sufficient to majorize the 
norms of the 2- forms Ri (the remainders) and 1-forms rji (incomplete ratios). In 
other words, one has to control only the division step, since multiplication by H is 
an operator whose norm is no greater than 

The division step is also rather transparent, its well-posedness being determined 
by the norm of the inverse Sylvester matrix -3^ from (7.2) and the norm of the non- 
principal terms \\H — L\\, where L is the principal homogeneous part of H. The 
group of affine transformations of the complex plane naturally acts on the space 
of all Hamiltonians of degree n + 1 not affecting the critical values of H. By an 
appropriate transformation of this group, one can always achieve the normalizing 
condition 113^^ H = 1 that is a condition on the principal homogeneous part L. 
The problem on bounding the norm \\A\\ is reduced therefore to studying how the 
magnitude of the non-principal coefficients of H — L may affect configuration of the 
critical values of H subject to the above normalizing condition on the principal part 
L. The problem can be explicitly solved for the univariate polynomials, implying 
an answer in the hyperelliptic case as well [NYOlb]. Moreover, one can show 
that if all critical values of a bivariate polynomial transversal to infinity coincide, 
then necessarily the polynomial H must coincide with its principal part, being thus 
homogeneous, eventually after a suitable parallel translation in the (a;, j/)-plane. 
The inequality between the non-homogeneity of a bivariate polynomial and the 
spread of its critical values is still unknown, see [NYOlb] for partial results. 

7.6. Preliminary conclusion. The tangential Hilbert problem for generic 
Hamiltonians, gets reduced to the question about the number of isolated zeros 
of linear combinations of functions satisfying together a system of linear ordinary 
differential equations with rational coefficients. 

The procedure of derivation of this system is very transparent. In particular, 
it can be written in the hypergeometric form (7.5) with explicit bounds on the 
norms of the corresponding matrices A and B. These bounds in turn imply that 
when reduced to the Fuchsian form (4.3), the system will have the residue matrices 
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Aj bounded (as usual, in the sense of the norms) in terms of the inverse distance 



When this inverse distance tends to infinity (which corresponds to confluence 
of singular points), the resulting bounds on the norms of the residues in 

(4.3) explode. However, this explosion is of a very specific nature: the monodromy 
group of the system with one or several confluent singularities, remains the same. 
In particular, the spectral data of the residues remain bounded. 

In the subsequent lectures we will find out how far away are these conditions 
from sufficient conditions allowing for an explicit solution of the tangential Hilbert 
problem. 

Lecture 3. Quasialgebraicity of function fields 

Starting from this moment, we will pursue the same path towards the tangential 
Hilbert problem, but this time in the opposite direction. Namely, we will establish 

conditions on Fuchsian systems guaranteeing that their solutions are similar to 
algebraic functions, in particular, admit explicit bounds for the number of isolated 
zeros. 



The main objects of study in this section are functional fields obtained by 
adjoining one or several analytic (in general, multivalued) functions to the field C 
of complex numbers (or slightly more generally to that of rational functions C(t)). 
Such fields admit filtration (grading) by degrees. 

The goal is to obtain conditions on the field (in terms of properties of the 
adjoined functions) guaranteeing that the question on the global number of isolated 
zeros of functions from this field can be algorithmically (effectively) solved. An 
accurate definition will be given at an appropriate moment, after explaining all 
technicalities pertinent to the problem. We begin by examples illustrating the 
goals. 

8.1. Algebraic functions. The field of rational functions C{t) in one variable 
t possesses the following obvious but nevertheless remarkable property: 

1. any element f{t) = p{f)/q{f), p.q E C[f.] from this field has a well defined 
degree deg / = max(deg j;. deg g) (assuming that the representation is irre- 
ducible), and 

2. the number of isolated zeros of / on the whole pro jective line CP^ = CU{oc.} 
is no greater than deg / (actually, equal to it if counted with multiplicities). 

In other words, there is a direct relationship between the combinatorial complexity 
of representation of / in the field (i.e., the number of field operations necessary 
to produce / from constants and the independent variable t), and its analytic 
complexity measured by the number of isolated zeros. 

This example can be easily generalized by considering fields generated by one 
or several algebraic functions. 

Let fi{t), i = 1, . . . ,n he algebraic functions of one variable, defined implicitly 
by the polynomial equations 
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with respect to a;,. Consider the ring C[/i, . . . , /„] formed by polynomial combi- 
nations of the functions /i, and the corresponding field of fractions C(/i, . . . , /„). 
Both consist of analytic multivalued functions ramified over a finite point set on 
the projective line CP^, though the number of distinct branches of every function 
is finite (no more than di ■ ■ ■ dn). 

One can define unambiguously the degree of functions in this ring and the 
respective field as the degree d of the minimal representation 

n 

Proposition 3. The total number of isolated zeros of a polynomial combina- 
tion f{t) = P{fi{t), fn{t)), f e C[/i, ...,/„] of degree ^ d on all branches of 
this function, does not exceed d ■ di ■ ■ ■ dn- 

Proof. This is an immediate corollary to the Bezout theorem applied to the 
system of algebraic equations 

Pi{t,xi) = 0, 

< 

_ P{xi,...,Xn) = 0. 

The same proof actually works for a more general case of functions fi defined by 
a system of algebraic equations Pi{t, xi,. . . , Xn) = 0, i = I, . . . ,n, involving all 
functions simultaneously. □ 

Remark. The function fo{t) — t is clearly algebraic, and if required, we can 
always assume it being among the collection of the functions fi, thus avoiding 
particular cases and awkward notation. This agreement will allow us to assume 
that all functional fields are extensions of the field C{t). 

8.2. Existential quasialgebraicity. Some parts of the above construction 
can be reproduced in a completely general context. Let U C C be an open domain 
(for simplicity assume it to be bounded, connected and simply connected). 

Consider an arbitrary collection of n functions F = {fi{t), . . . , fn{t)} analytic 
in U. They define the ring C[F] and the corresponding field C{F) of functions 
meromorphic in U. As before, for any function / € C[P] one can define its degree 
as the minimal possible degree of the polynomial P G C[xi, . . . ,Xn] realizing the 
given function f{t) = P{F{t)) = P(/i(t), . . . , /„(t)) (there can be algebraically 
dependent functions among the generators). This grading extends naturally for 
rational combinations from C{F). 

Proposition 4. For any compact K <i= U there exists a counting function 
C = Ck : N ^ N, taking only finite values Cxid) < +oo for any finite d, such that 
the number of isolated zeros of any function f in K can be at most Cxid): 

degf ^d ^ #{teK: f{t) = 0} < CK{d). 

Proof. First we show that under the assumptions of the Theorem, the number 

of isolated zeros of any linear combination fc(t) = ^ Cifi{t) with complex constant 
coefficients ci, . . . ,c„ G C, is bounded in any compact K uniformly over all such 
linear combinations. Indeed, without loss of generality we may assume that the 
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functions /j are linear independent — this does not affect the supply of all linear 

combinations. Next, it is sufficient to consider only combinations with coefficients 
on the unit sphere, satisfying the equality = 1- The functions fc for such 

c are all different from identical zero, hence each of them has only a finite number 
of isolated zeros in the compact K (accumulation of roots to the boimdary of K is 
forbidden). Now the standard semicontinuity arguments using compactness of the 
unit sphere, prove that the number of zeros of all fc is uniformly bounded. 

To deal with arbitrarily polynomial combinations, wc can treat them as linear 
combinations of monomials /"{t), \a\ ^ d, reducing the general case to the already 
studied one. □ 

Remark. One can easily recognize in this demonstration some minor varia- 
tions on the theme already exposed in §1. Of course. Proposition 4 follows from 
Theorem 2, since the parameters Ca can be considered varying over the compact 
sphere. Here we could explicitly avoid dealing with functions vanishing identically 
for some values of the parameters, simplifying considerably the proof. 

8.3. Comparison. Two above finiteness assertions, Propositions 3 and 4, dif- 
fer in two important instances: 

1. the bounds on roots of algebraic functions are global, i.e., valid on the max- 
imal domain of definition of the functions from the field C(i^), whereas the 
bounds on the roots of arbitrary analytic functions in general depend on the 
choice of the compact K and can in many cases blow up as the compact 
approaches the boundary of the maximal domain U ; 

2. the bounds on roots of algebraic functions are given by an explicit formula 
involving some basic parameters defining the field, whereas the function 
Cxid) is totally existential (see below). 

The nature of the counting function Cxid) from Proposition 4 remains totally 
non-effective. One can easily construct examples of functions (even entire functions) 
such that the growth of Cxid) will be arbitrarily fast [IY96]. 

8.4. In search of quasialgebraicity: reappearance of Picard Vessiot 
extensions. Our goal is to provide sufficient conditions on the functions fi guaran- 
teeing that the corresponding field will be similar to the field obtained by adjoining 
algebraic functions. This condition, still understood informally until made precise 
in §10, will be referred to as quasialgebraicity of the function field. The accurate 
definition is postponed since it involves some technical details. 

However, even prior to giving any accurate formulation, the class of function 
fields among which one could hope to find nontrivial cases of quasialgebraicity, can 
be substantially restricted. 

First, the generating functions /i, ■.-,/« must be multivalued (ramified). In- 
deed, a single-valued function having at most polar singularities on the projective 
line CP^ (recall that we are looking for global bounds, hence the functions should 
be defined globally), must be rational. Any field generated by such functions, is a 
subfield of C(t) and hence we get nothing new. 

On the other hand, if one of the functions fi has an essential singularity on 
CP^, then by classical theorems of complex analysis this function near such point 
must take infinitely many times almost all values, hence one can easily construct 
an individual polynomial combination having infinitely many roots accumulating to 
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the essential singularity. This precludes quasialgebraicity whatever exact meaning 
it may be assigned. 

Thus any field (or ring, what is almost the same for our purposes) exhibit- 
ing nontrivial quasialgebraicity, must consist of functions ramified over some finite 
set S C CP^. As above, globality means that functions should be analytically 
continuable along any path avoiding the ramification locus S. 

The possibility of analytic continuation along paths (and loops) introduces an 
additional structure, the monodromy group action. Choose arbitrarily a nonsin- 
gular point a ^ S. Then any element from the ring C[P] can be identified with 
the full analytic continuation of its germ at a. Denote as in §5 the monodromy 
operator associated with a loop 7 G 7ri(C \ 17, a) by A.^. Then it would be natural 
to assume that C[P] is closed (invariant) by analytic continuations, that is, 

A^/i GC[/i,...,/^] V7e7ri(CPKi:,a). (8.1) 

Moreover, since the filtration of the ring C[P] by degrees should be well-defined, it 
must be preserved by analytic continuations, that is, continuation of a polynomial 
/ = P(/i, . . . , /^) (of minimal degree among all polynomials representing the given 
/) along any loop should again be a polynomial of the same degree. In particular, 
the C-linear span of the germs fi, . . . , in the space of all analytic germs must he 
invariant by all monodromy operators. In other words, there should exist invertible 
square matrices My such that 

A.,(/i,...,/„) = (/!,...,/„) -M^. 

In the same way as in §5, this implies that the functions f\,---,fix must be 

solutions of a linear ordinary differential equation with single- valued coefficients on 
CP^ \ S. As before, allowing essential singularities of the coefficients would im- 
mediately relinquish any control over nonaccumulation of roots of solutions, hence 
we arrive to the following important conclusion: To he quasialgehraic, the field 
C(P) = C(/i,...,/^) must he a Picard-Vessiot extension of C or C{t) obtained 
by adjoining solutions of a linear ordinary differential equation with rational coef- 
ficients. 

This fact, in particular, implies that adding to the generating tuple of functions 
(/i, . . . , /„) their derivatives of orders ^ n — 1 will make the field C{t, /i, . . . , fn) 
a differential field (closed by differentiation). Without loss of generality we can 
assume that such completion was already done and we deal with the differential 
field C{X) (and the corresponding ring C[X]) obtained by adjoining to C all entries 
of a fundamental matrix solution X{t) for a system of first order linear ordinary 
differential equations with rational coefficient matrix, 

X{t) = A{t)X{t), A{t) e Mat^xM(C(t)). (8.2) 
In the particular case we are discussing, one can take X to be the Wronskian matrix 
of the collection {/J, Xij{t) = /j'~^^(t), i,j = l,...,n. 

8.5. Quasialgebraicity of Fuchsian systems: examples and counterex- 
cmiples. As in the case of single-valued functions, a precondition for quasialge- 
braicity is nonaccumulation of roots of solutions to singular points. Singularities of 
linear systems (8.2) can be easily described: they occur only at the poles of the co- 
efficients matrix A{t) or rather at the poles of the matrix-valued differential 1-form 
A{t)dt on the Riemann sphere CP^. The corresponding classical theory [Har82] 
distinguishes between two types of singularities: 
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• those exhibiting at most polynomial growth of entries of the fundamental 

matrix X{t) and its inverse X~^{t) and called regular singularities, and 

• those exhibiting abnormally fast (faster than polynomial) growth of solu- 
tions, called irregular singularities. 

Remark. In order to measure growth rate of multivalued functions near a 
ramification point, they should be restricted on a sector bounded by two rectilinear 
rays with the vertex at this point. Otherwise one can construct a curve approaching 
the singular point while spiraling around it in such a way that the growth in terms 
of the distance to the singularity will be arbitrarily fast even for the most innocent 
multivalued function Int. 

The dichotomy between regular and irregular singularities is closely related 
to dichotomy between poles and essential singularities for single- valued functions. 

Consider a neighborhood of a singular point t^ E assuming for simplicity that 
t^ = 0. Let M be the monodromy operator associated with a small loop around 
the origin, so that 

A^X(t) = X{t)M, det M 7^ 0. (8.3) 
Let A e Mat^X([*(C) be any matrix logarithm of M, so that 

exp 2mA = M. 

Then the multivalued matrix function t^ = exp{A\nt) also has the same mon- 
odromy, 

= cxp(AA^ Int) = exp(yl(lnt + 27rz)) = t^M = Mt^. 

Therefore the matrix ratio H{t) = X{t)t~'^ is single-valued in a small punctured 
neighborhood of 0. If the origin is a regular singularity, then H{t) has at most a 
pole and hence by choosing a different valuation of the logarithm and replacing 
by for sufficiently large natural k, one can make H{t) holomorphic at the 

origin. By explicit elementary formulas for matrix exponents, one can derive from 
this a local representation of entries of the matrix X (t) in the form 

k,\ 

with the exponents A ranging over the spectrum of A, the natural k being no 

greater than the maximal size of Jordan blocks of A and hk\{t) holomorphic at 
t = 0. Notice the remarkable coincidence with (4.1) and Exercise 1.2: the latter 
implies that at least in the situation when all functions hk\ and all eigenvalues of A 
are real, isolated zeros cannot accumulate to the regular singularity at the origin. 

If, on the other hand, the singularity is irregular, then at least some entries of 
the matrix H{t) must exhibit essential singularity at the origin. In the same way 
as with single-valued function, one can in this case construct functions from C{X) 
that would have infinitely many roots accumulating to the origin. 

Example 8.1. The linear system 

j ±1= t~'^Xi, 
\X2 = Q 



(8.4) 
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generates the field C{X) containing a function f{t) = exp(— 1/f) — 1 whose zeros at 
the points tm = (27rim)~^, m = ±1, ±2, . . . accumulate to two essentially singular 
points at t = and t = oo along the imaginary axis. A 

Thus occurrence of irregular singularities destroys any hope to achieve quasial- 
gebraicity, and we are left with the class of Picard-Vessiot extensions exhibiting 
only regular singularities. 

There is a simple sufficient condition guaranteeing that a singular point of the 
linear system (8.2) is regular. By the Fuchs theorem [H£ir82], if A{t) has a simple 
pole (of the first order) at a point f*, then this point is a regular singularity. Such 
singularities are called Fuchsian. 

The inverse to the Fuchs theorem is in general not true: there exist regular non- 
Fuchsian singularities. But a system exhibiting regular singularity at t^, by a m,ero- 
morphic (locally near i«) linear transformation X{t) R{t)X{t) can be reduced 
to the system having a simple pole at f*. Thus on the level of local meromorphic 
equivalence there is no difference between Fuchsian and regular singularities. 

For a globally defined system exhibiting several regular singularities at finitely 
many points S = {ti,. . . ,td,} on the projective line, one can ask whether there 
exists a globally meromorphic (hence rational) linear transformation simultaneously 
taking all regular singularities into Fuchsian ones. The problem (that constitutes a 
part of the so called Riemann-Hilbert, or 21st Hilbert problem) turns out to be very 
delicate, the result depending essentially on the structure of the monodromy group 
of the system, and not always admitting solution, as shown recently by A. Bolibruch 
[Bol95, Bol90]. However, from the classical result by Plemelj [Ple64, For91] it 
follows that if an additional singular point is allowed to be created anywhere, then 
the answer is always positive and a rational matrix function R{t) can be found such 
that Y{t) = R{t)X{t) satisfies a system of linear ordinary differential equations 
having only Fuchsian singularities on the whole projective line CP^. Clearly, if 
t G C(X) (which we may always assume without loss of generality, as noted above), 
then C{X) = C{Y) and hence when discussing quasialgebraicity, one can deal with 
Fuchsian systems only. 

A Fuchsian system with d finite singular points ti, . . . ,td can always be written 
in the form 

a;eC", e Mat„xn(C), t e C, ti,...,tdeC, 

explicitly indicating the corresponding residue matrices Aj (we return to the initial 
notation X for the dependent variables). Thus the natural problem arises. When 
the Picard-Vessiot extension C{X) constructed by adjoining all components of the 
fundamental matrix solution of a, Fuchsian system (8.5), is quasialgebraic? Note 
that this field depends not on the choice of the matrix solution X{t) but rather on its 
"logarithmic derivative" A{t) = X{t)X~^{t) which is a rational matrix function. 
The Fuchsian system (8.5) is determined by its dimension and the collection of 
algebraic data {Ai, ti, i = 1, . . . , d}. The bound for the number of zeros should be 
given in terms of these algebraic data. 

Example 8.2. The simplest class of Fuchsian systems is that having only two 
singularities (one simple pole of the matrix A{t) on the whole line CP^ is impossible 
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since the sum of all the residues, including the one at infinity, must be zero). By a 

conformal transformation of the independent variable the two points can be placed 
at t = and t = oo. The corresponding system will then take the Euler form, 

X = j-X, AeMat^x^(C). (8.6) 

The Euler system (8.6) can be immediately integrated: X{t) = t^. The asso- 
ciated Picard Vossiot extension has the form that is already familiar: 

C(JC) =C({t^ln*^"^OA,fc), A e Spec A, fc < m(A), (8.7) 

where m(A) > 1 is the maximal size of Jordan cells corresponding to the eigenvalue 
A. 

Roots of real functions from the field C(X) cannot accumulate to the origin, if 
all eigenvalues of A are real (Example 1.2). We will see that the "Euler field" for a 
system with only real spectrum is indeed quasialgebraic. 

On the other hand, the field C{t^^, t~^^) associated with the 2 x 2-Euler sys- 
tem with non-real spectrum {±^/^}, contains the function cos In t that obviously 
has infinite number of roots accumulating to f = 0. 

This suggests that quasialgebraicity of the fields C{X) should be somehow 
related to spectral properties of the residue matrices. -4 

8.6. Counting zeros of multivalued functions globally. Since functions 
from the Picard- Vessiot extensions are in general multivalued, counting their zeros 
"on the whole Riemann sphere" should be supplied with a precise meaning. 

Example 8.3. The function f{t) = — 1 for a real irrational A has infinitely 
many different roots tk — , q = exp /c G Z, lying on different branches of this 
function. 

In order to single out a branch, it is required to choose a simply connected 
domain C C \ i7. Usually this choice is done by slitting the complex plane along 
segments or rays with cndpoints in the singular locus S. In the above example, 
slitting C along the negative semiaxis yields a simply connected domain U with the 
following property: any branch of the function t^ — 1 in this domain has at most 
[A"""^] roots in it (the same bound holds also for —U). Here and below [aj stands 
for the integer part of a real number a. This explicit bound is global in the sense 
that U and —U together cover C \ Z'. < 

However, the choice of the simply connected domain U can affect even finiteness 
of bounds on the number of zeros. 

Example 8.4. If A is a non-real munber, then the roots fall on the logarithmic 
spiral (since \q\ ^ 1 in this case). If U is chosen by slitting C along positive or 
negative semiaxis, then the number of isolated roots on any branch will still be 
bounded by [(ReA)-iJ < -hoo for A ^ iM. 

On the other hand, one can choose a simply connected spiral domain containing 
the origin on the boundary, so that it would simultaneously contain an infinite 
number of the points tk = q'' ■ 

Though the spiral slits are not algebraic (even not analytic) curves, this example 
can be easily modified to construct polygonal simply connected domains containing 
as many roots of tlu; fmiction f'^ — 1 as necessary, provided that the number of sides 
of the polygon can be arbitrary. < 
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These examples motivate the following definition. 

Definition 6. An analytic multivalued function /: CP^ — > C ramified over a 

finite set E C CP^ is said to admit a global bound on the number of isolated zeros, 
if there exists a natural number N < +00 such that the number of isolated zeros 
of any branch of / in any open rectilinear triangle T c C \ S has no more than A'' 
isolated zeros, the bound being uniform over all such triangles. 

Remark. Having this property, one can place an explicit upper bound on the 
number of isolated roots of / in any simply connected domain U bounded by finitely 
many algebraic curves of known degrees. The number of curves, their degrees and 
the number of ramification points will explicitly enter the answer together with the 
number TV. 

Indeed, the domain C \ S can be triangulated into finitely many triangles as 
above. The number of simply connected components of any intersection UHT with 
each triangle of the triangulation can be easily majorized by the Bezout theorem, in 
terms of the above data. The total number of such components should be multiplied 
by N to construct the required upper bound. 

Example 8.5. The function t^ — 1 admits a global upper bound for the number 
of isolated roots for all A with Re A ^ 0. For A = •\/— T any real interval (0, a), a > 0, 
contains infinitely many roots of this function. -4 

Remark. In practice, however, we will always choose a convenient system of 

rectilinear slits as in the textbooks on complex variables: if, for example, an upper 
bound for the number of zeros is known in both [/ = C \ M.^ and —U = C\ 
then the number of zeros in any triangle not containing the origin, does not exceed 
the maximum of these two bounds, since such triangle cannot intersect both the 
positive and negative semiaxcs simultaneously. 

9. Digression on computability 

9.1. Hierarchy of computability. The discussion in the previous section 
addressed the issue of globality of bounds on the number of zeros of multivalued 
functions. Now we pass to a brief discussion of com,putahility of these bomids. This 
would lead us again to functions defined by differential equations with algebraic 
right hand side parts. 

The bounds for fimctions of the field C(/i, . . . , /^), besides being dependent on 
the degree d of these functions, should obviously depend on the field itself, that is, 
ultimately, on the choice of the generators /j. While computability of functions of 
natural arguments is a well developed area (several notions of computability can be 
distinguished, see below), dependence on functional parameters is definitely outside 
the scope of any theory (e.g., the number of roots of a function is not a computable 
characteristic per se). 

An intermediate position occupy "computable" functions of one or several real 
or complex variables. Here we take the most restrictive attitude, legalizing only 
polynomial functions, the absolute value | • | of a real or complex variable, and the 
integer part [-J of a real variable. 

Returning to functions of one or several natural arguments, one can classify 
them in an "increasing order of computability" into: 
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1. (general) recursive functions, whose values can be computed for any combi- 
nation of arguments in a finite time by some deterministic algorithm, Turing 
machine etc.; 

2. primitive recursive functions, that can be defined by one or several iterated 
inductive rules of the form 

(t>{k + l,m) =^{k,m,(j){k,m)), /c G N, m € N", (9.1) 

assuming that the functions and = ^(1, •) are already known; 

3. elementary functions, given by finite compositions of powers, exponents, 
algebraic operations etc. 

Example 9.1 (tower functions). Consider the function T{k,n) of two integer 
arguments, defined by the recursive rule 

r(0,fc)=fc, r(n + l,fc) = 2^^"''=). (9.2) 

These rules define tower functions (iterated exponents): r(n, •) is a tower of height 
n and for each particular n is an elementary function of the second argument k: 

r(l,fc) = 2^ T{2,k)=2^\ Ti3,k)=2^"\ ... 
However, as a function of the first argument, it is not elementary. M 

Example 9.2 (Ackermann generalized exponential) . Consider the function j4(^;, 
of three natural arguments, defined by the recursive rules and initial conditions as 
follows, 

A{z + l,x,y+l)= A{z, X, A{z + 1, x, y)), 

A{0, x,y)=y + l, A{1, x, 0) = x, A{2, x, 0) = 0, (9.3) 

A{z,x,0)=0, yz^2. 

These rules define each function (j) = A{z+1, •, • ) for any particular z unambiguously 
provided that the function = A{z, •, • ) is already defined, by the simple recursive 
rule (j){x, y + 1) = ^{x, <j){x, y)). In other words, the set of conditions defines A for 
all nonnegative combinations of arguments. 
One can easily check that 

A(l,x,y) = X + IA hl^x + w, 

A{^,x,y) = " 

A{2, x,y) = x + x + -- -+ x = xy, 
A(3, X, y) = XX ■ ■ ■ X = x^ , 

This suggests that as a function of the first argument, A is not primitive recursive. 
The reason is that, unlike in the simple recursive rule (9.1), the function A is applied 
to itself (this does not prove that A is not primitive recursive, but the fact remains 
true and can be rigorously demonstrated). It can be also shown that as a function 
of the first argument, the Ackermann generalized exponential grows faster than any 
primitive recursive function. M 

9.2. Transcendental functions defined by algebraic data. The above 
brief discussion suggests that in order to speak about computable bounds depend- 
ing on several analytic functions as parameters, these analytic functions must them- 
selves be defined in terms of finitely many integer, real or complex parameters. 
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The simplest example of such functions are polynomials (their coefficients play 

the role of the parameters) and, slightly more generally, algebraic functions. How- 
ever, this example does not allow to produce nontrivial quasialgebraic fields. 

Diff'erentiation of algebraic functions leaves them algebraic. On the contrary, in- 
tegration (taking primitives) in general results in transcendental functions. Another 
possibility of constructing transcendental functions from algebraic ones is exponen- 
tiation. However, both primitives and exponentials are only particular cases of 
functions defined by ordinary differential equations with algebraic coefficients (e.g., 
x{t) = exp J fit) dt is a solution to the differential equation x = f{t)x). This looks 
as a most general mechanism transforming algebraic input data into transcendental 
output functions. 

Thus we arrive to the following natural conclusion: in order to speak about 
quasialgebraicity of the function fields C(/i, . . . , /^), the generating functions fi{t) 
must be defined by polynomial ordinary differential equations or systems of siich 
equations. Then by computability of any bound we would always assume that 
this bound can be expressed in terms of the (real or complex) coefficients of these 
equations and/or integer parameters (dimensions, degrees, etc), so that as fimctions 
of real or complex parameters these bounds are polynomial, while being elementary 
(or at worst primitive recursive) functions of the remaining integer variables. 

Example 9.3 (principal). Suppose that the functions fi{t), i = to- 
gether satisfy a system of polynomial ordinary differential equations of the form 

Xi= ^ Ciax"", i = l,...,iJ., a = {ai,...,aij,) gZ^, 

\a\^m 

with complex coefficients Ci^ € C. Then the expressions 

Ci = (max|ci„|)^('''™), C2 = N{i,,m) ■ V |ci„| 

are both computable characteristics of the collection {/i, . . . , /^} provided that N 
is an elementary or at worst primitive recursive function of two integer arguments 
m (the degree) and /x (the dimension of the system). M 

9.3. Restricted computable bounds. Explicit computability does not im- 
ply (neither does it assume) globality of the bounds. In the same way as the 
parameters defining the functions fi may enter the answer, sometimes numeric 
characteristics of the domain (in our settings, always a triangle T C C \ 17) may 
enter into the expressions. 

Definition 7. Let 17 C C be a domain (usually polygonal) eventually contain- 
ing some of the singular points inside. We say that the field C{X) is quasialgebraic 
in U, if the number of zeros of any function from this field can be bounded in any 
triangle T C U \ E, uniformly on all such triangles, but the bound may depend on 
U. 

The most commonly occurring form of this dependence is through the distance 
between dU and the singular locus E. 

Remark. Since the independent variable ranges over the Riemann sphere CP^ 
and the point t = 00 may well belong to the singular locus, the distance from dU 
to S should be defined in such cases as minimum of the above distance to the 
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finite part of E and the number mit^gu \t ^\ measuring the "distance from dU to 
infinity" . 

10. Quasialgebraicity and uniform quasialgebraicity 

10.1. Quasialgebraicity of Picard— Vessiot extensions: accurate for- 
mulation of the problem. Let X(t) be a fundamental matrix solution of a Fuch- 
sian system of n linear ordinary differential equations (8.5), and C(X) (rcsp., C[X]) 
the field obtained by adjoining all entries of this matrix to the field C (rcsp., the 
ring of all polynomial combinations of these entries). As was already noted, we 
always assume that <C{X) contains the subfield of rational functions C(t). 

Definition 8. The field <C{X) is called quasialgebraic, if the number of isolated 
roots of any function / G C{X) of degree k in this field, in any triangle T free from 
singular points tj of the Fuchsian system (8.5), is bounded by a number depending 
only on: 

1. the degree k = degc(x) /; 

2. the dimension n and the number d of (finite) singular points (as an elemen- 
tary or at worst primitive recursive function); 

3. the entries of the residue matrices Ai, . . . ,A4 (in an algebraic way); 

4. on the coordinates of the singularities ti,. . . ,td (also in an algebraic way). 

In fact, in all cases when quasialgebraicity of the field will be established, the 
bounds would depend on the complex (matrix) parameters in a very simple way, 
via the residual norm of the rational matrix function A{t), 

i?(A(.))= max (10.1) 

j = l,...,d 

In a similar manner, dependence of the bounds on the position of the singularities 
tj will be expressed via the (inverse) spread of these points, 

p{E) = mBD^{\ti - tj\-^ , \ti\}. (10.2) 

(this number is large only when some of the singular points approach each other). 

10.2. Uniform quasialgebraicity. According to the above definition, the 
field is quasialgebraic if the simple algebraic data (parameters) defining it can be 
used to produce an explicit bound on the number of zeros of any function from this 
field. 

Some of these parameters are always relevant. There is no question why the 
dimension and the degree must necessarily enter any bound on zeros. It may require 
some eff'orts to see that other parameters can also aff'ect the answer. 

Example 10.1. The Euler system 

xi=at~^xi, X2 = 0, a € N, 

defines a field containing the function f{t) — t°- — 1 that is of degree one in this 
field (though, as a polynomial, it has degree a). The number of isolated zeros of 
the function / in a triangle can be as large as [|J . < 

This example shows that the magnitude of eigenvalues of the residue matrices 
clearly affects the number of zeros. 

It is not very diSicult (see the subsequent sections) to construct upper bounds 
(global or not) that would involve rational expressions of the parameters (in this 
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case, entries of the residue matrices) that have poles, making the bounds exploding 
for certain combinations of parameters. On the other hand, there are no visi- 
ble reasons for appearance of the infinite number of zeros for these values of the 
parameters. It requires considerable efforts to show that in terms of the residue 
matrices, the bounds can be given by polynomial expressions, which is equivalent 
to expressing them in terms of the residual norm (10.1) as above. Finally, it is 
the most difficult part to prove that under certain additional but rather natural 
assumptions, the bounds can be given uniformly over all configurations of singular 
points. 

10.3. Notes, remarks. We conclude this highly informal discussion by sev- 
eral remarks, also of a very general nature. 

As stated in §4, the ultimate goal of the theory is to establish a constructive 
bound for the number of zeros of Abelian integrals (tangential Hilbert problem). 
These integrals were shown to belong to certain Picard-Vessiot extensions. Appli- 
cation of the methods explained in subsequent sections requires repeated algebraic 
manipulations with the integrals, therefore bringing into play the whole field of 
functions. Moreover, on the final stage of the construction not one but rather 
several Picard Vcssiot fields arc considered simultancoiisly, while carrying out in- 
duction in the number of ramification points. This explains why the existence of 
upper bounds on the number of zeros (quasialgebraicity) is defined as a property of 
the corresponding functional fields (rather than elements constituting these fields). 

There was also a special reason for choosing the definition of quasialgebraic- 
ity without attempting to speciiy explicitly the "counting function" measuring the 
number of isolated zeros. This was done primarily because the bounds that one can 
obtain on this way are enormously excessive. Several nested inductive construc- 
tions immediately produce tower-like bounds even from the modest exponential 
contributions on each inductive step. 

Finally, we would like to note that the Fuchsian representation is not the only 
possible. Actually, the Picard-Fuchs system of differential equations for Abelian 
integrals is written in the hypergeometric form 

{tE + A)X{t) = BX{t) 

determined by two constant matrices A,B & Mat„xra(C), and it would be natural 
to require that quasialgebraicity were expressed in terms of the norms ||^||, 
rather than in terms of the respective residues after transforming to the Fuchsian 
form. This type of bounds is not proved (if it were, this would imply constructive 
solution of the tangential Hilbert problem, as explained in §7.3). 

Lecture 4. Quantitative theory of differential equations 

11. Bounded meandering principle: explicit bounds on zeros of 
functions defined by ordinary differential equations 

The preceding section contained motivations for introducing the notion of 
quasialgebraicity. In this section we survey some results of constructive (though not 
always global) nature, bounding the number of isolated zeros of functions defined 
by differential equations. 
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11.1. Linear n order equations with bounded analytic coefficients. 

The basis for all other considerations is a classical theorem by de la Vallcc Poussin 
concerning solutions of linear ordinary differential equations with bounded coeffi- 
cients, not necessarily polynomial or even analytic. It gives a sufficient condition 
for a linear ordinary differential equation of order n with real bounded coefficients 
guaranteeing absence of solutions with more than n — 1 isolated roots on a given 
real interval 7 c M. Such equations are called disconjugate on I. In order to stress 
the difference with the complex case, we denote the independent variable by s, and 
consider a linear equation 



+ ai(s) + • • • + a„_2(s) y" + a„_i(s)y' + a„(s)2/ = (f f .1) 



on the real interval I = [sq, Si] of length r with real bounded coefficients: |a;c(s)| < 
Cfe < oo for all s G I. 



then any -smooth function /(s) satisfying a, linear equation (11.1) may have at 
most n — 1 isolated roots on I , counted with multiplicities. 

This result can be seen as a generalization of the Sturm nonoscillation theorem 
for equations of order greater than 2. The inequality can be slightly improved, see 
[Lev69]. 

This (simple) statement implies a number of corollaries. First, subdividing 
any interval into sufficiently short segments satisfying (11.2), and adding together 
the bounds for each segment, one can obtain an explicit bound on zeros valid on 
any real interval where the coefficients a,(s) are explicitly bounded. Next, one can 
consider equations of the same form (11.1) with complex-valued coefficients and 
solutions. There is no sense to count zeros of such solutions (since in general there 
are no roots), instead an interesting question is to estimate their topological index, 
the variation of argument Arg/(s) between the endpoints of the real interval. By a 
simple modification of the method used in the proof of Lemma 8 one can obtain an 
upper bound for this topological index, again in terms of the magnitude (uniform 
upper bounds for the absolute value) of coefficients of the equation (11.1), see 
[Yak99], as follows. 

Lemma 9. Variation of argument of any solution f{t) of the equation (11.1) 
with complex-valued coefficients all hounded by C in the absolute value, along the 
segment I of length \I\ = si — sq is no greater than 



If the equation has analytic coefficients Oj G C'(t^) defined in some domain 
U C C, then the number of complex isolated zeros inside any polygonal domain 
can be majorized using the above described bounds for the topological index and 
the classical argument principle. As a result, the following inequality for the number 
of isolated zeros can be obtained (we return to the initial notation t for the complex 
independent variable). 



Lemma 8 ([dlVP29]). If 




(11.2) 



7r(n + l)(l + 3C|/|). 
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Theorem 7 ([Yak99]). Suppose fi{t),.. ., fn{t) solve 

+ ai(t)j/("-i) + • • • + an-i{t)y' + a„(i)y = 0, t€U, 

ai&0{U), \ai{t)\ ^R, i = l,...,n. 

Then for any triangle T C U of perirneter £ any linear combination f 
with complex constant coefficients Ci, has no more than 

|(n+l)(l+«) 

isolated zeros in T . 

Remark. Analyticity of solutions of the equation (11.3) follows from analytic- 
ity of its coefficients. However, the theorem apparently can be modified to cover also 
equations with meromorphic coefficients, provided that the solutions are analytic 
in U (the so called apparent singularities). 

11.2. Systems of linear equations: Novikov's counterexample. Moti- 
vated by potential applications to Abelian integrals and quasialgebraic fields, we 
would like to find generalizations of Lemma 8 or Theorem 7 for systems of linear 
equations, in hope to majorize the number of isolated roots of any linear combina- 
tion of coordinates cixi{s) H + c„a;„(s) for a system 

x{s) = A{s)x{s), se/cM, Vse7 ||A(s)|| < ii, 

in terms of ^ = |/| and R. The following counterexample due to D. Novikov 
[NovOlb], shows that it is impossible, neither for real nor for complex systems. 

Example 11.1. Let ti,...,td be collection of different numbers (real or com- 
plex) from a real interval / or a triangle T (consider only the last case for simplicity). 
Denote a{f) = e{t — ti)...{t — td) and let £ be so small that |a(t)| + |a(f) + a^(t)| < 1 
for any t€T. 

Solution (f)i = exp( J a{t)dt) of the linear differential equation of the first order 
Xi = a{t)xi, has no zeroes at all. However, the derivative (1)2 = (f>i = a{t)(j)\ has 
the same zeroes as a{t) and satisfies the equation ^2 = (a + a^)^i- 

Together the pair {(j)i,(j)2) satisfies the linear 2 x 2-system 

d 

xi=axi, X2 = {a + a^)xi, a{t) = s'^^{t — tj), (H-S) 

1 

whose coefficients are bounded by 1 everywhere in T by the choice of e. However, 
the second component has the specified number d of isolated zeroes there, where d 
can be arbitrarily large. < 

This example suggests that when dealing with systems of linear equations, 
instead of general analytic (variable) bounded coefficients, only rational or even 
polynomial coefficients were allowed, so that 

d 

A{t) = ^kt^^ e Mat„xn(C), (11.6) 

fe=0 

and the magnitude of the matrix coefficients, e.g., the total norm R = Ylik W-^kW 
used to produce upper bounds. 



(11.3) 
(11.4) 
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11.3. Reduction of a system to an equation: obstructions. The above 
example is not very surprising since the standard procedure of reducing a linear 
n X n-system 

x = A{t)x, X = {xi,...,x„) gC^ (11.7) 

to one nth order linear equation is discontinuous with respect to parameters. In- 
deed, without loss of generality one can assume that it is the number of zeros of 
the first component t/ = Xi in a domain U C C that interests us. Differentiating xi 
by virtue of the system, we see that the derivatives are linear combinations: 

y^''\t)=qf,{t)-x{t), k=l,...,n,..., (11.8) 

where the analytic (co)vector functions Qkit) = {qk,i{t), ■ ■ ■ ,qk,n{t)) € C"* are 
determined by the recurrent rule 

qo(t) = (l,0,...,0), q„+,{t) = q^it) + q,,{t) ■ Ait) (11.9) 

and do not depend on the choice of the trajectory x{t). (We use the right multipli- 
cation to stress that q^it) are row vector functions). Consider the field k = M{U) 
of functions meromorphic in U and the linear n-space k" over this field. 

The functions qj.(-) are vectors in this space. Consider the linear subspaces 

{0} C io C ii C • • • C Ln-i C L„ C • • • C k", (11.10) 

each Lk being spanned by qo,qi, ■ ■ ■ ,qk over k. If all inclusions are strict, the 
dimensions of L/. over k must strictly increase, while being bounded by n = dimif k". 
Hence the ascending chain of subspaces (11.10) must stabilize no later than after 
n steps, i.e., for some jj, ^ n the inclusion must be non-strict, G L^-i, which 
implies a linear identity over k = A4{U), 

q^ + ai{t)q^_i + ---+a^-i{t)q^+a^{t)qo = 0, ai€k = M{U). (11.11) 

This implies that y = q^ ■ x satisfies the equation (11.3). 

Though the recurrent formulas generating the vectors 9i, ^2, • • • are explicit, it 
is impossible to control the magnitude of the coefficients aj(-), whatever that may 
mean, from this construction. This can be seen in a simpler settings when the field 

k is replaced by C. 

Example 11.2. Consider a sequence of vectors qo,qi, . . . , in the n-dimensional 
space over M or C, that grows at most exponentially in the sense of the norm: 

ii^fcii s^c^ c>o. 

For the same dimensionality reasons as above, no later that at the nth step a 
linear dependence must occur, allowing to express some qk as a linear combination 
of the preceding vectors go, • • • , qk-i, 

-qk + aiqk-i-\ \-akqo = 0, e C, (11-12) 

similarly to (11.11). However, coefficients of this dependence arc out of control and 
can be arbitrarily large. To see this, consider the situation when the angle between 
qk-i and Lk-2 is very small but nonzero so that Lk-2 £ Lk-i, while qk belongs 
to Lk-i and is orthogonal to Lfe_2 in La-i- Even assuming some lower bounds on 
ll^fcll will not help to improve the situation in this case. 

The situation becomes completely different under any of the following two 
additional assumptions. 
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Example 11.3. Assume that all the vectors qk actually belong to a lattice Z" c 

C" and their norms arc bounded from above as in the preceding example. Then the 
coefHcients of the linear combination (11.11) can be explicitly bounded from above: 
indeed, if qu is a linear combination ofqo,. . qk-i, then the coefficients a, in (11.12) 
can be found by solving a system of linear nonhomogeneous algebraic equations with 
integer matrix of coefficients and integer free terms. By the Cramer rule, solutions of 
this system can be obtained (after elimination of redundant equations and assuming 
that solutions indeed exist) as ratios of appropriate minors (determinants of some 
square submatrices) . These minors are integer numbers, explicitly bounded from 
above by virtue of the assumptions on the norms \\qi\\. Hence in each ratio the 
numerator is bounded from above, whereas the denominator is no smaller than 1 
in the absolute value. This clearly implies an upper bound on all the coefficients 

Example 11.4. Assume that the vectors qk € C" are obtained by iterations of 

a linear map P: C" — > C". Then the bounds on the norms will be automatically 
satisfied with c = ||-P||. The first linear combination between the vectors still 
can have very large coefHcients, exactly as in Example 11.2. However, instead 
of looking for the first combination, we may continue until the step number n. 
Since any operator P is a matrix root of its characteristic polynomial x{z) = z" + 

a\z"~^ + h ttn-iz + a„, the linear combination of the form (11.12) in this case 

can be obtained via the identity x{P)'lQ = 0- Coefficients of the characteristic 
polynomial admit an upper bound in terms of the eigenvalues of P, each of which 
is no greater than the norm ||P||. Finally, we conclude that a linear combination 
(11.11) in this case (iterations of a linear map P in the finite-dimensional space 
C")) can be constructed with coefficients Ui explicitly bounded in terms of 
and n. M 

11.4. Reduction of a system to an equation: commutative algebra 
versus linear algebra. The operator q{t) q{t) + q{t) ■ A{t) from the n-dimen- 
sional C(i)-linear space into itself, describing the iterations (11.9), is not linear 
over the field C{t), since the derivative d/dt is not C(<)-linear operator. On the 
other hand, considered over the field C, the space of polynomial vector-functions is 
not finite-dimensional unless deg A{t) = (in which case the degrees are uniformly 
bounded). Thus neither of the above methods can work. 

However, a solution can be found in terms of the commutative algebra. It allows 
to treat in a similar way both linear and nonlinear systems. Before proceeding 
further, we give a useful technical definition. 

Definition 9. The norm, of a polynomial p € C[xi, . . . ,Xn] is the sum of 
absolute values of its coefficients: 



One can easily verify that in addition to the usual triangle inequality, this norm is 
multiplicative, \\pq\\ ^ \\p\\ \\q\\ for any two polynomials. The norm of a derivative 
\\dip\\ (in any variable) can be easily bounded, \\dip\\ < degp \\p\\. 

Generalization of the construction from §11.3 is very simple. 

Example 11.5 (basic). Consider the linear system (11.7) with a polynomial 
matrix of coefficients A{t) as in (11.6). The rule (11.9) defines a sequence of 
polynomials qk{t,x) = qk{t) ■ x G C[t,x], that are always linear in the variables 




(11.13) 
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X = (xi, . . . , x„): the recursive rules allow to estimate the degrees degj qk{t, x) and, 
if necessary, their norms. 

Instead of the linear subspaces L}~ over the field k = <C{t) in this case, consider 
the polynomial ideals Ik = {qo, ■ ■ ■ ,qk) C C[f, a;]. The ascending chain 

{0} C 7o C 7i C . . . C 7„ C . . . C C[t, x] (11-14) 

of these ideals must eventually stabilize in the sense that some inclusion Ii-i C Ic 
becomes non-strict (an equality). This follows from the fundamental fact that the 
polynomial ring C[t,x] is Noetherian and any ascending polynomial chain in it 
eventually stabilizes. 

The above stabilization condition means that ±qc € h-i, hence for appropriate 
polynomial coefficients hi, . . . ,he G C[t, x] 

e 

qi + ^hiqe-i = 0. (11.15) 

i=l 

A priori, the polynomial coefficients hi can depend on a; in a nontrivial way. 
However, since all polynomials cjk arc linear in x (homogeneous), one can truncate 
the identity (11.15) retaining only constant terms of hi (of degree in x) and 
construct a new identity of exactly the same form (11.15) but with coefficients 
ai{t) = hi{t,0) from the univariate polynomial ring C[t\, 
i 

qi + '^aiqi-i = 0, ai € C[t], i = 1, . . . ,£. (11.16) 

i=l 

The identity (11.16) obviously means that the function qo{t,x) = x\ satisfies 
the linear ordinary differential equation 

2/(^)+^ai(t)t/(^-*) =0, ai{t)=hi{t,Q)&C[t\, i=l,...,L (11.17) 

i=l 

In order to apply the results form §11.1, one has to estimate the absolute values 
|ai(t)| from above. In any disk of known radius {\t\ < R} this can be done if the 
decomposition (11.15) is explicitly known, in particular, if the following parameters 
are explicitly bounded from above: 

• the length £ of the chain (equal to the order of the resulting differential 
equation) ; 

• the degrees of the coefficients hi in t; 

• the norms i = 1, . . . , ^. 

•4 

Example 11.6. The construction involving chains of ideals, is not very degree- 
specific. Consider a system of polynomial ordinary differential equations, 

Xi = Pi{t,x), i = l,...,n, PieC[t,xi,...,Xn]- (11.18) 

Let / be a polynomial combination f{t) = Q{t,x\{t), . . . ,Xn{t)) for some Q = 
Q{t, x) € C[t, x]. Consider the infinite sequence of polynomials qo, qi, . . . , q^, • • • € 
C[t,a;] formed by iterations of the Lie derivative, i.e., the recursive rule 

^ dqk , -sr^ dqk „ /n 1n^ 

qo = Q, (lk+i = -Q^+2^Q^,-Pi- (11-19) 

i=l 
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Let Ik = {qo,qi, ■ ■ ■ ,qk) C C[f, x] be the polynomial ideals generated by the 
first k + 1 polynomials from this sequence in the polynomial ring C[i,a;]. They 
obviously form an ascending chain, Ik C Ik+i- Since the ring C[t,a;] is Noetherian, 
this chain must stabilize and hence for some natural £, 

qe + hiqe-i + ■ ■ ■ + hi-iqi + htqo = 0, hi€C[t,x]. (11.20) 

As in the linear case, this identity implies a polynomial relationship between the 
unknown function y = Q{t,x(t)) and its derivatives up to order i. Unlike the 
linear case, this time the coefficients hi may depend on x explicitly, so that in 
addition to the data described in the previous example, one needs upper bounds 
for \xi{t)\ in [/, which is a non-algebraic piece of information. However, in many 
cases this information can be easily achieved (or even a priori known) and the 
problem reduces to getting the same information as in the linear case, namely: the 
length £ of the ascending chain of polynomial ideals, the degrees and the norms of 
the polynomials hi appearing in the representation (11.20). < 

11.5. Generalizations and improvements. It is not clear from the very 

beginning, how replacing ascending chains of linear subspaces by ascending chains 
of polynomial ideals may resolve the problems related to unboundedness of the 
coefficients of the decompositions (11.11) and (11.20) respectively, see Exercise 11.2. 
We explain it now. 

The first advantage of the suggested approach allows to treat systems depending 
polynomially on additional parameters, without the risk of producing bounds that 
blow up for certain values of the parameters. Actually, the difference between 
parameters and phase variables disappears almost completely. 

Example 11.7. Assume that the linear system (8.2) with a polynomial matrix 
A{t) (11.6) depends on additional parameters A = (Ai, . . . , Ap) is a polynomial way: 
A e Mat„xn(C[f, A]). Then one can consider the chains of ideals in the bigger ring 
C[t, A] and only minor notation changes are necessary to construct a linear ordinary 
differential equation (11.17) whose coefficients will be in fact polynomial in t and A 
(note that the resulting differential equation is always monic: its leading coefficient 
before the principal derivative is 1). 

This polynomiality eliminates the danger that for some values of the parameters 
the coefficients of the derived equation (11.17) will blow up (which was earlier the 
case). < 

Another advantage appears as a generalization of Example 11.3. 

Example 11.8. Assume that the coefficients matrix A{t) of the system (11.7) 
is polynomial as in (11.6), and in addition all matrix coefficients Ai- have only 
integer entries: Ak € Mat„xn(^)- Then, since the polynomial qo(t,x) = xi also 
belongs to the subring Z[t,x], all subsequent polynomials qk will also have integer 
coefficients, and their degrees are growing not faster than linearly in k. The growth 
of the norms \\qk\\ can also be easily controlled. 

Suppose that the length i of the ascending chain of the corresponding ideals 
(11.14) is already known (in the univariate case it is relatively simple, see [NY99b]). 
Then one can explicitly compute an upper bound r for the degrees of polynomial 
coefficients hi in (11.15) in terms of i, the degree d = degA{t) and n = dimx. 
To find the polynomials hi themselves (or rather the univariate polynomials a, = 
hi{-,0)), it is possible now to use the method of indeterminate coefficients: writing 
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each Oj as J2]=o '^"io^^ substituting them into the identity (11.16), we obtain a 
system of non-homogeneous linear algebraic equations for the unknown variables 
{cij : i = 1, . . . , j = 0, . . . , r}, with integral coefficients matrix and the free terms 
column, all explicitly bounded from below. For the same reason as in Example 11.3, 
in this "lattice" case explicit bounds for ||aj|| can be immediately produced in terms 
of n, d and r. < 

The lattice polynomial system may look artificial, but in fact their appearance 
is natural. The explanation is given in the following principal example. 

Example 11.9 (universal system). All entries of the matrix coefficients Ak of 

the polynomial matrix function (11.6) can be considered as parameters and denoted 
by Ai, j = 1, . . . , (d + l)n^. With respect to these variables, the "universal matrix 
polynomial" Yl'i^o ^kt^ is in fact in the "lattice" Mat„xn(Z[i, A]), of known degree 
d+ 1: moreover, all coefficients of A{t) over the ring Z[t, A] are zero-one matrices. 

In a similar way one can start iterations from the "general linear form" go(ij x) = 
Y^=i f^iXi which becomes a polynomial of degree n + 1 with zero-one coefficients 
in C[a;, /3i, . . . , /3„] and add the string of the coefficients to the parameter list. 
This would allow to treat simultaneously isolated zeros of all nontrivial linear com- 
binations. M 

11.6. Lengths of ascending chains. Example 11.6 in fact proves that for 

any dimension n and degree d of the matrix polynomial A(t) from (11.6), one can 
majorize explicitly the number of isolated zeros of any linear combination ^ PiXi {t) 
in terms of n,d, and R — J2o ll^fcll provided that the length of ascending chain 
of polynomial ideals (11.14) is explicitly known. We claim that this length is a 
"computable" function. 

Notice that for each combination of n, we have a uniquely defined parameter 
space {A} — c(''+-'^)" (coefficients of the matrix polynomial A{t) and the initial 
linear form qo), and the uniquely defined chain of ideals in the polynomial ring 
C[t, X, A] generated by polynomials qk that in fact belong to the lattice Z[t, x, A]. 

The construction of the generators qk is absolutely explicit. The problem of 
verifying whether the next polynomial qk belongs to the ideal Ik-i generated by 
the previous polynomials, is constructive: there exists an algorithmic procedure 
allowing to get a positive or negative answer in a finite number of steps for each 
k. Thus the length £ is a "computable" (in the weakest sense) function of n and d: 
the above description can be transformed into an algorithm computing the value 
£{n, d) for any given combination of n and d. 

Coupled with the preceding discussion, this computability means that the prob- 
lem on the number of zeros of functions defined by polynomial ODE's, is algorithmi- 
cally solvable. However, the complexity of this algorithm turns out to be very high 
in the sense that the bound on the growth rate of the function £(n, d) is tremendous 
and higher than any elementary function. In the next section we address this ques- 
tion in details and show what other modifications are necessary in order to produce 
a "theoretically feasible" upper bound. 

12. Lengths of chains 

12.1. Descending chains of algebraic varieties. The problem on ascend- 
ing chains of polynomial ideals, belongs entirely to the realm of commutative alge- 
bra. However, there is a parallel geometric problem that admits a considerably more 
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transparent solution. Moreover, the construction used in the "geometric" proof, can 

be adjusted (after introducing appropriate technical tools) to the algebraic case. 
Consider a strictly decreasing chain of complex algebraic varieties, 

= XqD XiD X2D ■■■ D XkD ■■■ (12.1) 

where each variety Xk is given by a finite number of polynomial equations in C" 
of degree no greater than dk (no restrictions are imposed on the number or the 
structure of these equations). Without loss of generality we may assume that the 
sequence c?i, c?2, • • • is non-decreasing. 

By the same Noetherian property, the chain (12.1) must terminate after finitely 
many steps at some variety X^. The problem is to compute the length £ from the 
dimension n and knowing the bounds for the degrees dk- 

If the degrees dk are all bounded, then the problem belongs to linear algebra, 
since the polynomial equations defining the varieties Xk will in fact constitute a 
finite-dimensional linear space, and its dimension will be the natural bound for 
the length of the chain (12.1). In order to avoid technical troubles when talking 
about computable dependence on infinite input data {di,d2, ■■■,}, we will consider 
finite-parameter examples. The most important are the linear growth case, when 

dk = d + k, fc = l,2,... (12.2) 

or the exponential growth case 

dk=d^, A; = 1,2,.... (12.3) 

In both cases the natural number d is a parameter, and in both cases t should be 

majorized in terms of n and d. 

12.2. Lexicographically decreasing sequences of words. The advantage 
of the "geometric" problem on chains of complex algebraic varieties, is that it 
can be reduced to a purely combinatorial problem on lexicographically decreasing 
sequences of words. 

Recall that any algebraic variety can be uniquely represented as the union of 

irreducible algebraic subvarictics of different dimensions varying from (isolated 
points) to n — 1 (hypersurfaces). Thus with any X C C" one can associate a 
vector v{X) e Z" with n integer nonnegative coordinates , . . . , u^,v^), where 

J/* — v^{X) stands for the number of irreducible z-dimcnsional components of X. 
Denote by -< the lexicographic order on letting 

(ry"-S . . . , J.0) -< {D"-\...,D°) 

if and only if for some k = 1, . . . ,n 

^n-l ^ ^n-l^ _ _ _ ^ = but I/^-^ < d''-\ 

The following elementary observation is crucial. 
Lemma 10. IfXCYc C", then u{X) -< i/{Y). 
Proof. This holds since: 

1. any irreducible component of X should belong to an irreducible component 
of Y, and 

2. if A c B is a pair of irreducible varieties, then dim74 < dimS and in the 
case of equal dimensions necessarily A = B. 
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In other words, when passing from y to X each irreducible component either com- 
pletely survives, or is split into a number of other components of lower dimen- 
sions. □ 

The fact that a descending chain of algebraic varieties must stabilize, follows 
now from the following purely combinatorial claim. 

Proposition 5. A lexicographically strictly decreasing chain 

iyiyu2> ^ z/fe > , i^k&^+, (12.4) 

must be finite. 

Proof. For n = 1 the claim is obvious, since any decreasing sequence of 
nonnegative integers must be finite. For an arbitrary n, the first "letters" z^^"^ € 
Z+, k — 1,2, . . . , form a non- increasing sequence (it must not necessarily be strictly 
decreasing). However, no more than a finite number of values is taken. Along 
any interval of constancy of the first letter, the tails (i/^~^, . . . , f°) also form a 
lexicographically strictly decreasing sequence in Hence the length of each 

such segment is finite by the induction assumption, and the length of the whole 
chain is finite as the sum of finitely many finite numbers. □ 

This proof, being extremely simple, can be supplied by quantitative estimates 
of the lengths under additional assumptions on the norms of the words. Denote 
= H h the norm on Z^. 

Example 12.1. Assume that the norms of the words U). forming the chain 
(12.4) are bounded: 

\\i^k\\ ^ d + k, yk = 1,2, . . . {d, a natural parameter). (12-5) 

Then the length of the chain as a function of n and d is bounded by a general 
recursive (but not primitive recursive) function. The explanation is as follows. 

Let f{n, d, i) be the maximal length of the chain under the additional constraint 
that the first "letter" of the first word is no greater than i, ^ i- Then the 

restricted function /(n+1, •, i + 1) can be expressed via /(n, •, •) and /(n + 1, • • • 

Indeed, the length of the initial segment of the chain, on which the first letter 
maintains its initial value i, can be at most f{n,d,d), since the tails start with a 
word whose first letter can be at most d. After the first segment is exhausted, the 
remaining part of the chain begins with a word of length n, whose first letter is at 
most i — 1, and with the restriction on the norms of the words as follows, 

\Wk\\ ^d+{k + f{n,d,d)), 

if the words of the remaining chain are numbered after the drop of the first letter. 
This is tantamount to replacing d by f{n,d,d), therefore we obtain the recurrent 
inequality 

fin + 1, d, i) ^ f{n, d, d) + f{n+l,d + f{n, d,d),i- 1). (12.6) 

Coupled with the boundary conditions 

fin,d,0)^fin-l,d,d), f{l,d,i)^i, 

this determines the upper bound for / completely. Notice the remarkable similarity 
between the recurrent rule (12.6) for / and that for the Ackermann generalized 
exponential (9.3). < 
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The growth rate determined by the above recursive function, is enormous: in 
terms of n, the dimension of the ambient space, the function / grows faster than 
any primitive function. On the other hand, this bound is essentiahy sharp: for lex- 
icographically ordered chains it is rather obvious, for chains of polynomial ideals it 
was proved by G. Moreno [MS92, MS91], for chains of algebraic varieties the claim 
is apparently also true though not written anywhere. In any case the bounds based 
on such estimate of the length of chains should not be considered as constructive, 
even theoretically. 

However, under rather mild additional assumptions the bounds can be improved 
very considerably, in fact to become elementary functions of n, d. 

12.3. Dynamically generated chains. The chain of ideals (11.14) is gener- 
ated by a dynamical system. More precisely, the rule (11.19) for the general case of 
a polynomial vector field (11.18) is a Lie derivative of the ring C[i, x\, and the ideals 
Ik are generated by iterated Lie derivatives of the initial polynomial qo G C[t,x]. 

A parallel construction for chains of algebraic varieties also involves a dynamical 
ingredient [NY97, NY99b]. Namely, instead of the general chains (12.1) with 
the only restriction on the degrees of the varieties Xk, one should consider chains 
generated by a discrete time dynamical system in C". 

Let F : C" C" be a polynomial map of some known degree d and X C C" an 
algebraic subvariety. Consider the decreasing chain {X^} defined by the recursive 
rule 

Xo = X, Xk+i=XnF-'{Xk), A; = 0,1,2,.... (12.7) 
In other words, 

Xk = xn F-^{x) n F-^{x) n • • • n f-''{x), 

where F'' stands for the k times iterated map F o ■ ■ ■ o F and F~''{-) = (F'^)~^(-) 
denotes the corresponding preimage. 

Dynamically Xk can be described as the set of points a G C" such that the 
point a itself, together with its F-orbit F{a), . . . ,F''{a) of length k, belong to X. 
This immediately implies that 

F{Xk X Xk+i) C Xk-i \Xk, fc = 1, 2, . . . (12.8) 

(the differences occurring above consist of initial conditions of orbits that jump off 
the variety X exactly after the specified number of steps). 

The condition (12.8) means that the difference Xk \ Xk+i cannot be too large 
compared with the preceding difference Xk-i \ Xk- For an arbitrary polynomial 
map this is not true, but under the additional assumption that F preserves di- 
mensions of semialgebraic sets (i.e., takes curves into curves and not into points, 
though eventually creating singularities, and the same in higher dimensions), one 
can conclude that the dimensions of the differences Xk \ Xk+i are non-increasing: 

dim(Xfc\Xfe+i) <dim(Xfe_i \Xfe), fc=l,2,.... (12.9) 

12.4. Dynamically generated chains stabilize fast. Consider a strictly 
descending chain of varieties (12.1) satisfying the additional condition (12.9) on 

the dimensions. Then the associated words I'k = i^{Xk) E Z", in addition to the 
lexicographic decrease (12.4), display stronger monotonicity properties. 

Consider the sequence of "heads" (as opposed to "tails"), obtained by trun- 
cating the words Uk to their first s symbols, [uk]s = ('^fe"^)" " " j^k~^) ^ 
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all k = 1,2, For any lexicographically non- increasing sequence, the sequence 

of heads of any fixed length s will be again non-inereasing (by definition of the 
lexicographic order). Yet this monotonicity can be non-strict in general. 

However, under the additional assumption (12.9) any truncated sequence of 
heads {[i^fe]s}^i for the sequence Vk = z^(^fe) will be strictly decreasing in the 
following sense: 

Vs = l,...,n-1 ^ Hs > y [J^k-ih = Wk]s = ■ ■ ■ (12.10) 

In other words, as soon as the first equality between truncated words occurs, the 
rest of the chain will have the same heads. Indeed, [ffe_i]s = [i'k]s if and only if 
all irreducible components of Xk-i and Xk of dimensions n — s and higher, are the 
same, which means that the difference Xk-i \ Xk is at most {n—s— l)-dimensional. 

Such type of descent ensures much faster convergence. Indeed, the first letter 
stabilizes after no more than steps, after which the problem is reduced to 

that for the tails, which are words of length n — 1. The corresponding inductive 
inequality is very simple. 

Example 12.2. Denote hy g{n, d) the maximal length of decreasing sequence of 
words of length n growing no faster than linear in the sense of the norm, W^kW ^ d+ 
k, with d being an integer parameter as in (12.2), under the additional assumption 
(12.10). Then the above observation implies that the length until stabilization of 
the first letter is at most d, whereas the sequence of tails of length n — 1 starts from 
the word of norm no greater than d + d = 2d and hence its length by the induction 
assumption is no greater than g{n— 1, 2d), so that finally 

g{n,d) ^d + g{n-l,2d), g{l,d)^d. (12.11) 

This gives an upper bound for the length, 

g{n, d) (2" - l)d, 

which is much better than in the general case considered in Example 12.1. A 

Example 12.3. In a similar way the exponential bound (12.3) can be treated. 
In this case the inductive inequality analogous to (12.11) takes the form 

g{n,d)^d + g{n-l,d'^), (12.12) 

which gives g(n, d) as a tower function of height n. Though not elementary, this is 
obviously a primitive recursive function of both arguments. -4 

These simple examples illustrate the algorithmic complexity of the problem on 
lengths of lexicographically descending chains. For descending chains of algebraic 
varieties the answer follows in fact from Example 12.3, since the degrees of equations 
defining the algebraic varieties Xk grow exponentially: degF*^ = d^ , where d = 
deg-F (without loss of generality we may assume that degX — d with the same d). 
The bounds for the norms ||!/(Xfe)|| follow from a version of Bezout theorem due to 
J. Heintz [Hei83]. 

The most technically difficult case is that of ascending chains of polynomial 
ideals, mainly since there is no uniqueness in the primary decomposition of such 
ideals, hence one cannot associate a word v{I) with an ideal /, counting the number 
of primary components of various dimensions. However, the components of the 
maximal dimensions are correctly defined and may be counted, which allows to 
implement a similar inductive proof. Details can be found in [NY99b] . 
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13. Restricted qucisialgebraicity of Picard— Vessiot fields 

Quasialgebraicity of a function field was defined as a property allowing for 

counting zeros globally. However, for technical reasons we need a weaker notion 
of restricted quasialgebraicity. Assuming t/ C C being a simply connected (usually 
polygonal or circular) domain containing no singular points on its boundary (but 
eventually some singularities inside U), we can restrict functions from C{X) on 
U. The result, denoted by Cu{X), is another functional field, again consisting of 
multivalued functions. 

Definition 10. We say that C{X) is restricted quasialgebraic in U, if the up- 
per bounds on the number of isolated zeros in any triangle T C U \ S can be given 
in the same terms as before, plus eventually some geometric parameters describ- 
ing the relative position of S and U, most often the distance between S and the 
boundary dU. 

13.1. Bounds in the disk for polynomial systems. The above discussion 

allows to analyze completely the polynomial case. 

Theorem 8. // the coefficients matrix A{t) of a linear system (8.2) is polyno- 
mial as in (11.6) with bounded matrix coefficients of known total norm R, then the 
corresponding Picard-Vessiot extension C{X) will be quasialgebraic after restriction 
on any disk C C of radius r. 

In other words, the bounds on the number of zeros of functions from C£)^(X) 
will be explicit but depending on r. As r ^ +oo, the bounds explode. 

Explosion of the bounds occurs for a very simple reason: the singular point 
at infinity is in general an irregular (and certainly non-Fuchsian) singularity which 
may be an accumulation point for isolated roots of polynomial combinations (see 
§8.4). 

13.2. Fuchsian system: reduction to the polynomial case. To apply 
the previous technique to a Fuchsian system (8.5), it is sufficient to introduce the 
new independent (complex time) variable. 

The matrix function A{t) from (8.5) can be reduced to the common denomina- 
tor, 

d-l d 

A{t) = ^^p{t), p=Y, p,t\ x{t) = n(* - 

i=0 i=l 

The corresponding linear system is orbitally equivalent to the polynomial system 

r.=^^*^^' •=^, (13.1) 

1 t = x{t), dr' 



* dz 



where r e C is the new time variable. The map 

t ^ T{t) = I 

Jo X{z) 

is defined on the universal covering surface of C \ i7 and takes explicitly bounded 
values away from E. The inverse map r i— > t{T) covers the complement C \ 1^ so 
that for any triangle T (s C \ i7 on distance e > from S one can find p such that 
T is covered by the image of the disk {|t| < p} C C. 
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Application of Theorem 8 to the system (13.1) impUes the following corollary 
(after some preliminary work). 

Corollary 3. The field C{X) constructed from the Fuchsian system, (8.5), is 
quasialgebraic in any domain containing no singular points from S. The bounds 
would depend on r = r{U) = dist(9U', E) = dist{U, S) and explode as r ^ 0+. 

13.3. Rational systems with apparent singularities. Actually, not all 

singular points arc dangerous. As was already remarked, Theorem 7 is not about 
the number of zeros of solutions inside a triangle T, but rather about the variation 
of argument of these solutions along the boundary of T. 

This means that the above theorems on restricted quasialgebraicity of Fuchsian 
systems away from their poles, can be reformulated as assertions on computability 
of upper bounds for the variation of argument along polygonal paths distant from 
the singular locus E. 

If the domain U has only apparent singularities inside, then all solutions of the 
system (8.2) are meromorphic in U and the order of their poles can be explicitly 
bounded in terms of the corresponding residue norms By the argument 

principle, this together with the bounds on the index along the boundary, implies 
restricted quasialgebraicity of Cu{X). 

This remark shows that it is the multivalucdness of solutions rather than any 
other circumstance, that is an obstruction for the global quasialgebraicity. In the 
next sections we show how multivalued functions can be treated, first locally and 
then globally. 

Lecture 5. Isomonodromic reduction principle and Riemann— Hilbert 
problem 

14. Isomonodromic fields: local theory 

14.1. Euler field: an example. The Eulcr field (8.7), obtained by adjoining 
to C all entries of the (multivalued) matrix function t^ = exp(Alni) solving the 
system (8.6), already occurred as a first example when quasialgebraicity could be 
expected. 

The necessary condition for quasialgebraicity is given in terms of the spectrum 
of the constant matrix A (the only residue of the system). We have already seen 
that if some of the eigenvalues are non-real, then infinitely many real zeros may 
easily accumulate to the origin. It can be relatively easily seen that under this 
assumption on the spectrum of A, such accumulation is impossible not only for 
real, but also for complex zeros. Actually, a stronger assertion holds: the number 
of isolated roots can be explicitly majorized. 

The following result (together with explicit bounds that in this case are rather 
accurate) can be found in [K Y96] . Consider a finite set of exponents S' C C of di- 
ameter diamS' = max;^^A'eS whose points may have non-trivial multiplicities 
v{\) e N, so that by definition '}2\es ^i^) — 

Theorem 9. If S CR, then the number of isolated roots of any finite sum 

^CfcAi^ln'^"'*, Xg S, k^u{X), CkxeC, (14.1) 

X,k 

in any triangle T C C \ {0} is no greater than #5 —1-1-2 diamS'. □ 
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This theorem impUes quasialgebraicity of the Euler field, since any element / 
from C[t^] can be represented in the form (14.1) with explicit control over diamS 
and #5 expressed in terms of \\A\\, dim A and deg/. 

The assumption of this theorem can be formulated in terms of the spectrum of 
the (only) monodromy operator. 

Corollary 4. If the spectrum of the monodromy operator M corresponding to 
a small loop around the origin, belongs to the unit circle {i.e., all eigenvalues have 
modulus 1), then the Euler field C{t^) is globally quasialgebraic, and the number of 
isolated zeros can be bounded in terms of the norm \\A\\. 

14.2. Restricted quasialgebraicity near a Fuchsian point. Consider a 
Fuchsian system (8.5) with a singular point = at the origin and all residues 
of norm ^ R. Let p be a (sufficiently small) positive number such that all other 

singularities are at least 2/3-distant and at most 1 / p-distant from the origin. 

Let M be a monodromy matrix corresponding to a small positively oriented 
loop around the origin (this matrix is defined up to a conjugacy) and 5 C C the 
spectrum of M (i.e., the collection of eigenvalues which is independent of anjrthing 

but the system and the singular point). 

Theorem 10. If S C {|A| = 1}, then the field C{X) is quasialgebraic in the 
disk Dp = {\t\ < p}. The bound on the number of isolated zeros can be given in 
terms of R and 1/p. 

The outline of the proof of this theorem occupies the rest of this section. A 
similar result was proved in [RY96] for linear equations of order n near a Fuchsian 
singular point. 

14.3. Joint fields. Consider a Fuchsian singular point ti = of the system 

(8.5) with the corresponding residue matrix Ai = B. Then the fundamental matrix 
solution X{t) in any disk D containing no other singularities, can be represented 
as 

X{t) = Y{t)t'^, teD, (14.2) 

where F is a meromorphic (single- valued) function oft in D. Replacing if necessary 
B hy B + rE with an appropriate r € Z, one can always consider Y{t) as a holo- 
morphic function at t = 0. Expanding all matrix products, we conclude that the 
field C{X) belongs to a bigger field C{Z,Y) spanned jointly by entries of the two 
matrix functions, Z{t) = t^ and Y{t). Whereas the field C{Z) is a Picard-Vessiot 
extension for the Euler system, this is not immediately obvious concerning the ex- 
tension field C(y). However, we claim that the field C{Y) is a subfield of a bigger 
Picard-Vessiot field for some other Fuchsian system having t = as an apparent 
singularity. Indeed, 

Y = Xt-^ - t-'^Xt-^ B = A{t)Y - t-^YB. (14.3) 

This is not a linear system for y as a matrix function to satisfy (since it involves 

the matrix multiplication from both sides). But if all entries of Y are arranged as 
one column vector, then (14.3) becomes a system of n'^ linear ordinary differential 
equations with rational coefficients, exhibiting a Fuchsian singularity at t = (ac- 
tually, at all other poles of A{t) as well). The residues of this larger-size system 
can be explicitly constructed and their norms bounded from above. 
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The joint field C{Z, Y) possesses a property that was already used for "single" 
extension fields. 

Lemma 11. Variation of argument of any element from, the joint field C{X. Y) 
along a polygonal path distant from the union of singular loci, is a computable 
function. The same is true also for sufficiently small circular arcs around each 
singular point. 

Note that f = by constniction is an apparent singularity for this system. 
Thus the Picard-Vessiot field C(F) is quasialgebraic restricted on D, as follows 
from §13.3. 

The problem that we face now, is to prove the restricted quasialgebraieity in 
D of the joint extension field Cd{Z,Y) containing Cd(X) by virtue of (14.2), 
having already established restricted quasialgebraieity of each of the fields Cd{Z) 
and Cd{Y) separately and the fact that Y has only one apparent singularity in 
D. This will be done by reducing the question on (restricted) quasialgebraieity of 
Cd{Z,Y) to that for C{Z), using the fact that Y has a trivial monodromy in D. 
This reduction can be extracted from the papers by Petrov [Pet88, Pet90]. 

First we need some real analysis. 

14.4. Real closedness. Real part of an analytic function / on an open do- 
main U is not analytic (unless it is a constant). Yet the real part of the restriction 

on an analytic curve j G U can be extended to a neighborhood of 7 and some- 
times to the whole of U as an analytic function (of course, taking non-real values 
outside of 7). 

If J7 = D is a disk centered at the origin and 7 the real line M, then for any 
function / meromorphic in D one can take 

Re/(i) = i(/(i)+7(I)), Im/(f) = i(/(0-7(l)), (14.4) 

and then Re/ and Im/ will again be meromorphic in D and equal to the real 
(resp., imaginary) part of / on M. 

This observation allows to assume without loss of generality that all generators 
of the field C{Y) are real on the real axis. Otherwise, one should take their real and 
imaginary parts Re Y and Im Y as above (both are meromorphic matrix functions 
real on Mni?) and consider the field C(ReF, Imy) obviously containing C{Y). All 
properties of the field C{Y) are inherited by C(Rey. ImY). 

For a multivalued function / ramified over the origin t = G I?, taking its real 
or imaginary part on the whole real axis is an ambiguous operation because of the 
ramification: the real part on the positive semiaxis M+ can cease to be real after 
continuation on the negative semiaxis Yet one can often guarantee that the 
field as a whole is closed by taking real/ imaginary parts along any segment of the 
real axis. 

Lemma 12. Let A e Matnxn(R) be a real constant matrix and C{Z), Z{t) = 
t^ , the corresponding Euler field. Then for any function f G C(Z) its real or imag- 
inary part can be extended from the positive semiaxis M+ or the negative semiaxis 
M_ to functions g± again belonging to C{Z): 

V/ e C{Z) 3g± eC{Z): Re /|m± = 5± |m± , 

and the same for the imaginary parts Im/|K_,_ . 
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We will denote the functions g± = Re /kj. and their imaginary counterparts by 
Re± / and Im± / respectively. 

Proof. The field <C{Z) is independent of the choice of the fundamental so- 
lution, hence for each semiaxis ]R± one can choose a solution Z± of the Euler 
system tZ = AZ that is real on that semiaxis (recall that the matrix A is real so 
A{t) = t~^A is real-valued on R). Then it remains only to define g± using the 
identities 

Re(^c«ZS) =^(Rec„)Zg eC[Z], c„ G C, 

a 

and similarly for the imaginary part (with obvious modifications for the field of 
fractions C(Z)). □ 

Remark. The fact that the matrix A is real, is not a restriction: one can 
always consider a bigger Euler field of dimension 2n with the real block-diagonal 
matrix diag(Re A,lm.A). 

14.5. Variation of argument and zeros of imaginary part. The following 

elementary statement will play the central role in the constructions below. 

Lemma 13. ///: M D [a, 6] ^ C is a complex-valued function having no zeros 
on the real interval \a,b], then variation of argument of f along this interval is no 
greater than 7r(#{Im/ = 0} + 1). 

Note that here no extension from [a, b] is required, hence Im / stands for the 
usual imaginary part of the restriction. 

Proof. If | Arg/(fi) — Aigf{t2)\ > tt, then the imaginary part Im/ must 
vanish somewhere between ti and t2 on [a, b] . □ 

Actually, this result is true for any linear combination a Re / + /3 Im /. The 

function / may be allowed to have isolated zeros (but certainly not vanishing iden- 
tically), provided that zeros of Im / are counted with multiplicities. 

14.6. Isomonodromic reduction: the local case. Now everything is ready 
to prove that quasialgebraicity of C{Z) implies the restricted quasialgebraicity of 

C{Z,Y) in D = {\t\ < 1}, provided that Y{t) is single- valued in D and satisfies a 
Fuchsian system with bounded residues and all singularities e-distant from dD. 
Any function / e C{Z, Y) can be written in D as 

N 

/(i) = E^»Wy»(*)' ZiGC{Z), ViGCniY), i=l,...,N, (14.5) 
1 

where j/j are monomials from the field C{Y) and hence real on R, and the number 
of terms N as well as the degrees deg Zi , deg j/j (with respect to the corresponding 
fields C{Z) and C_d(F)) are exphcitly computable in terms of deg / with respect 
to C{X) c C{Z,Y). Indeed, for each generator Xij of degree 1 in C{X) such 
representation can be chosen with no more than n terms, each being bilinear in Y 
and Z. 

Assume that the field C{Z) is quasialgebraic. Then division by zi changes the 
number of zeros by no more than some known number. Hence when counting the 
number of zeros of / in any triangle T C D \ {0}, one can assume that zi{t) = 1, 
that is, f = yi+ ^iVi- 
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Also without loss of generality we may assume that T does not intersect the 
real axis (this can be also always achieved by subdividing it into smaller triangles 
if necessary), and belongs to the upper half-plane. Then the number of zeros of 
/ in T is no greater than the sum of four terms, variation of argument of / along 
the sufficiently small arc {\t\ ~ So <^ 1, Imt > 0} around the origin, and similar 
contributions from the semicircle {\t\ = 1, Imt > 0}, and two rectilinear intervals 
[—1, —So) C M_ and (eq, 1] C R+ respectively. 

The contributions from the two circular arcs arc computable functions by 
Lemma 11. As for the two rectilinear segments, the contribution of each of them 
is bounded by the number of zeros of Im / on them. By the assumption on the 
matrix function Y, Imj/j = 0, hence 

JV N 

Im± / = Im 2/1 + ^ ?/i Im = ^ 2/j Im± Zi . 

2 2 

By Lemma 12, the field C{Z) is closed by taking imaginary parts, hence for all i the 
functions Im-t Zi arc again in C{Z). Thus the question on the number of zeros of / 
represented by 14.5 in T is rcduc;cd to that for the number of zeros of two functions 
Im_|_ /, Im_ / € C{Z,Y) (one for each semiaxis), each of them involving less terms 
(at most N — 1). This allows to continue the process inductively, reducing the 
problem for zeros of / to that for some 2^ functions from <C{Z). On the last step 
the number of isolated zeros of a product ZnVn is majorized by the sum of zeros of 
each term (known since quasialgebraicity of C{Z) and C(F) is already established). 

Remark. The above described construction reducing quasialgebraicity of the 

joint field €.d{Z,Y) to that of Cd{X) if Y has trivial monodromy in D, is very 
similar to the standard differentiation-division scheme (based on the Rolle lemma) 
used to obtain bounds on the number of isolated zeros of real functions. Here the 
role of the differentiation is played by the operators Im±, whereas single- valued 
functions real on K play the role of constants killed by differentiation. 

14.7. Non-uniform quasialgebraicity. From Theorem 10 one can immedi- 
ately derive non- uniform quasialgebraicity of Fuchsian systems. Consider such a 

system with d singular points. 

Notice that the residue matrices Aj are invariant by conformal automorphisms 
of the independent variable, hence the residual norm of the corresponding matrix 
function R = maxi^i^....^ \ \Ai\\ is also invariant. 

Using such conformal automorphisms, one can always place any three poles of 
A{t) at any three points, say, 0, 1 and oo, but starting from the fourth pole, one has 
a nontrivial parameter characterizing the spread of singular points on the sphere. 
Let phe a small positive number such that: 

1. \U - tj\ > 2/9, i,j = l,...,d,i^ j; 

2. \U\ <l//9, i = l,...,d. 

Consider the monodromy matrices Mj corresponding to small loops going 
around tj . As before, their spectra Sj are uniquely defined. 

Theorem 11. If all spectra Sj belong to the unit circle, then the field C{X) is 
quasialgebraic. The bound for the number of zeros can be given in terms of R and 
P- 

Proof. Draw disjoint circles of radius p around each singularity and of ra- 
dius 1/p around the origin (this circle bounds a neighborhood of infinity on CP^). 
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Restricted on each circle, the field is quasialgebraic by Theorem 10. On the com- 
plement there arc no singularities, so after triangulation of this multiply connected 
domain one can apply Corollary 3. □ 

15. Uniform quasialgebraicity of Fuchsian systems 

In this section we briefly explain additional work to be done in order to obtain 
the bounds for quasialgebraicity, that would be independent on the relative position 
of singular points. In other words, we look for bounds that would remain explicit 
and uniform over p as the latter tends to zero. 

15.1. Isomonodromic reduction: the general case. The constructions of 
the previous section can be easily modified for a more general situation. Assume 
that two Fuchsian systems with fundamental solutions X[t) and Z{t) arc isom,on- 
odromic relative to a domain U, that is, they have the same singular locus S in U 
and the monodromy matrices Mj for all loops entirely belonging to U, are the same 
for the two. Next, assume that the residues (all of them, including those at singular 
points outside U) are all explicitly bounded. Finally, assume that U is & polygonal 
domain (say, a triangle) and the boundary dU is away from all singularities. 

Then one can claim that the two fields C(X) and C{Z) when restricted on U 
are both quasialgebraic or not quasialgebraic simultaneously. 

The proof in the case when all singularities in U fall on one straight line, is very 
similar to the local case. Namely, consider the matrix fraction Y{t) = X{t)Z^^{t) 
possessing trivial monodromy in U and embed C(X) into the joint field C{Z,Y). 
Writing elements of this joint field as Yl Ih^i ^nd applying the above algorithm of 
alternating division and taking the imaginary parts, the question on the number 
of zeros of / G C{X) can be reduced to that for several auxiliary functions from 
C{Z). 

15.2. Inductive strategy. The isomonodromic reduction principle as de- 
scribed above, would allow for an inductive proof of the uniform quasialgebraicity 
of Fuchsian systems if one could always construct a Fuchsian system that would be 
isomonodromic to a given one in a specified simply connected domain, while having 
no other singularities outside this domain. 

The inductive proof may look as follows. For Fuchsian system with only two 
singularities (Euler systems) the quasialgebraicity is known. Assume that it is 
already established for all systems with less than d finite singularities, and consider 
a system with d finite singular points forming the locus S C C. 

As was already noticed, it is impossible to make a conformal transformation 
placing all d points of S well apart from each other. However, one can always 
achieve a situation when all finite singular points form a set of diameter exactly 1 
inside the disk of radius 1 centered at the origin. In this case one can draw a line 
that is at least l/2(i-distant from all points of E and such that to each side of this 
line lies at least one (hence at most d — 1) point(s) of S. 

Now one can easily construct two polygonal domains each containing no more 
than d — 1 points of E, together covering the whole of E and with boundaries 
distant from both E and the infinity. 

Assume that for each such domain U a Fuchsian system can be found so that 
it will be isomonodromic with the given one in U, while still having the residual 
norm bounded in terms of the residual norm of the initial system. 
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Then application of the isomonodromic reduction principle would allow to re- 
duce the question on quasialgebraicity of the initial system in U to that for a 
Fuchsian system with ^ d — 1 finite singular points. By the inductive assumption, 
the latter question can be explicitly answered. 

The only assumption to monitor along this inductive process, is that on eigen- 
values of the monodromy operators. Clearly, one should assume that the spectral 
condition (on unit absolute values of eigenvalues) holds for all small loops around 
singular points. However, it is not sufficient, since this condition docs not survive 
the above surgery (cutting out part of the singularities inside U and pasting out 
the rest). Indeed, after replacing all singularities outside U by one singular point 
at infinity, we create a point whose local monodromy coincides with that of the 
boundary dU. Thus one must additionally assume that the monodromy along the 
boundary of U must also satisfy the spectral condition. 

Unfortunately, there is no way to predict how the partition into distant "clus- 
ters" of singular points will proceed when carried out inductively. Instead it is 
sufficient to assume that the spectral condition is satisfied for any simple loops 
(geometrically non-sclfintersccting closed Jordan curves). 

This conditional construction would prove the following theorem. 

Theorem 12. // the monodromy operators along all simple loops have only 
eigenvalues of modulus 1, then the Fuchsian system is uniformly quasialgebraic. 
Upper bounds for the number of isolated zeros can be given in terms of the number 
of singular points and the residual norm of the system, uniformly over all configu- 
rations of the singularities. 

However, the proof above gives only a general idea of how the actual proof is 
organized. The difficulties are of two kinds, technical and fundamental. 

An example of the technical problem is the isomonodromic reduction principle: 
it was formulated (and is actually proved) for a particular case when all singular 
points lie on just one straight line (which can then be identified with the real axis). 
However, it is sufficient for the purposes of the proof after suitable preparation of 
the initial Fuchsian system. 

The fundamental problem concerns the (im) possibility of constructing a Fuch- 
sian system with the prescribed monodromy. It was recently discovered that there 
exist obstructions to solvability of this problem, that have to be somehow circum- 
vented. In addition, one has to redress the proof of the corresponding positive 
results so that they would become constructive, yielding bounds for the residual 
norm of the constructed systems. 

16. Quantitative Riemann— Hilbert problem of matrix factorization 

16.1. Riemann— Hilbert problem: background. Given d distinct points 
ti, . . . ,td on the Riemann sphere CP^ and d invertible matrices Mi, . . . , Md satisfy- 
ing the identity Mi • • • M^ = E (the identity matrix), construct a Fuchsian system 
having singular points at ti, ... ,td and only there, for which the matrices Mj would 
be monodromy factors for some fundamental solution. 

This is the strongest form of the problem known as Hilbert 21st problem or 
the Riemann-Hilbert problem. Here some of many known results concerning its 
solvability. 

Plemelj theorem: For any collection of points and any monodromy matrices, 
one can construct a linear system with all but one singularities Fuchsian; the 
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last singular point is regular and can be made Fuchsian if the corresponding 
monodromy matrix is diagonalizablc. As a corollary, one can always con- 
struct a Fuchsian system with one extra singular point that would be an 
apparent singularity [Ple64, AI88]. 
Bolibruch— Kostov theorem: If the monodromy group generated by the ma- 
trices Mj is irreducible, then the Riemann-Hilbert problem is solvable [Bol92, 
Kos92]. 

Bolibruch counterexamples: There exists a reducible group generated by 
three 4 x 4-matrices, which cannot be realized as a monodromy group of a 
Fuchsian system [Bol90]. 

What is missing in all these formulations, is the possibility of placing any upper 
bound on the norms of the residue matrices of the corresponding Fuchsian system 
(provided that the problem is solvable). 

16.2. Surgery. In this section we formulate an analytic problem that is suf- 
ficient to solve in order to construct a system isomonodromic to the given one in a 
domain with all singularities away from the boundary. 

Recall that we have a Fuchsian system that already possesses the poles inside 
U at specified points and the monodromy around these poles is as required. What 
is necessary to do is to remove all finite singularities from C\U, leaving only one 
singular point at infinity. 

For simplicity we assume that U is the unit disk centered at the origin, and the 
annulus K = < \t\ < 2} is free from singular points (i.e., SdK = 0). Denote by 
X{t) the (multivalued) fundamental matrix solution. In general, X is multivalued 
even in K, however, for an appropriate constant matrix B the product W{t) = 
Xt~^ is a single- valued matrix function that is holomorphic and holomorphically 
invertible in K. 

Suppose that the matrix function W — Xt~^ , holomorphic and invertible in 
K can be represented as the matrix ratio of two other matrix functions, Ho{t) and 
Hoo{t), so that: 

1. Hq is holomorphic and holomorphically invertible in the disk Dq = {\t\ < 2}; 

2. Hao{t) is holomorphic and holomorphically invertible in the disk D^o = 
{\t\ > 5}, including the point t = 00; 

3. on the intersection K = DqC) , 

X{t)t-'' = Ho\t)H^{t), teDoHD^. (16.1) 

Then the two expressions, Xo(t) = Ho{t)X{t) defined in Dq and Xoo{t) = 
Hoo {t) defined on D^o , agree on the intersection and hence their "logarithmic 
derivatives" B(,{t) = X(,{t)X^^{t) and B^{t) = X^{t)X^^{t), coincide on K and 
together define a single- valued meromorphic function B{t) on the entire sphere 

The poles of B{t) can be easily described: due to the holomorphic invertibility 

of i?o,oo(i) they can occur either at the poles of A{t) that are inside Do, or at 
t = 00. In both cases the residues can be easily computed: for any tj € S Dq, 

B{t)=H^{t^).-^.H^\t^) + ... , 
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where the dots stand for terms holomorphic at tj. In a similar way, at infinity we 
have 

B{t) = ifoo(oo) • J ■ H-\oo) + 0{l/t^), 

where 0(l/t^) stands for a holomorphic matrix function of the corresponding 
growth. This means that in fact 

^W=E^' Bj = Ho{tj)AjH^\tj), (16.2) 

and Bj = ffoo(oo) B if-i(oo). 

The Fuchsian system with the matrix of coefficients B{t) given by (16.2), would 
serve our purposes, since its solution HqX has the same monodromy in Dq as the 
initial matrix solution X. 

In order to complete the proof, one would have to obtain upper bounds on the 
matrices Hj = Hci{tj) and their inverses for all singular points tj & S Do- 
This would imply an upper bound on the residual norm of the matrix function B(t) 
and finally would allow for the inductive proof as explained in §15.2. 

16.3. Birkhoff— Grothendieck factorization. Unfortunately, finding factor- 
ization (16.1) satisfying all properties above, is impossible. 

Example 16.1. One obstruction to holomorphic factorization can be imme- 
diately seen. Consider the determinant w{t) = det W{t): this is a holomorphic 
invertible function in the annulus K, and variation of argument of this function 
along, say, the middle circle of the annulus is an integer number i>. If this num- 
ber is nonzero, then either detHo or detH^ must have zeros and/or poles in the 
respective domains Dq resp., D^o. 

If n = 1, i.e., all matrices are of size 1x1, then one can always achieve 
factorization of the form 

W{t)=HQ^{t)t''H^it) (16.3) 

with 1 X 1-matrices Ho,Hoo holomorphically invertible in the respective domains. 

•4 

Both the positive and the negative assertions present in the above example, 
admit generalization for the general n-dimcnsional matrix case. 

Theorem 13 (see [GK58, GK60]). A matrix function W {t) holomorphic and 
invertible in the annulus K, can he factorized as follows, 

W{t) = t'^ H^{t), G = diag(i/i, . . . , i/„), Ui G Z, (16.4) 

with the matrix factors Hq,Hoo holomorphic and invertible in Dq resp., Doo- 

The integer numbers Vi, called partial indices, are uniquely determined {the 
same for all representations with the above properties) . 

The decomposition (16.4) can be used for surgery, if we incorporate the term 
t'^ into Hoo{t)- Then the matrix function Y{t) defined as HqX in Dq and as 
t'^Haot^ in Doo, satisfies a linear system of ordinary differential equations with 
only Fuchsian singular points in Dq and a regular non- Fuchsian point t = oo. The 
principal Laurent part of the matrix B{t) = YY~^ at y = oo can be easily described: 
the multiplicity of the pole at infinity is bounded in terms of v = \\G\\ = maxj Ivj], 
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and the magnitude (norm) of the Laurent coefficients of B at f = oo is bounded in 
terms of the Taylor coefficients of order ^ of the matrix Hco and its inverse. 

Actually, the form of the matrix t'^H^it) is not important: in order to have 
a system with a regular singularity at infinity, it is sufficient to have a matrix 
factorization 

1/ 

W{t) = Ho^H^, H^{t)^^ = C±,o{t) + C±,jt^, (16.5) 

j=i 

with holomorphic invertible matrix Hq and meromorphic invertible matrix function 
H^{t) having (together with its inverse H^) the only pole at t = oo of order < v 
with the Laurent (matrix) coefficients C±,j. 

Remark. The assertion of Theorem 13 is certainly not the strongest known. 
The matrix function W{t) can be defined only on the middle circle {\t\ = 1} of 
the annulus K and be rather weakly regular on it, still the factorization will be 
possible then, with the terms iJo, holomorphic invertible inside (resp., outside) 
the circle. Some non-circular contours can be also allowed. 

16.4. Modified surgery. One can easily modify the structure of the induc- 
tive construction above to cover the case of systems having only Fuchsian finite 

singular points and a regular eventually non-Fuchsian singularity at t = oo. Such 
system has the matrix of coefficients that can be written always as 

j = l ^ i=0 

and the residual norm for such systems should be defined as 

P(-)ii = E 11^^-11 + Eii^^ii- (16-7) 

3 i 

The surgery described in §16.2, using the factorization (16.5), allows to pass 
from one system from such class to another system from the same class, having 
no finite singularities outside Dq. In order to carry out inductively the bounds 
for zeros, one has to majorize the magnitude of all Laurent coefficients of the new 
system in terms of the norm R — ||A(-)|| of the initial system. 

16.5. Bounds. Suppose that a system from the class (16.6) of explicitly bound- 
ed norm R = ||^(-)|| has no singularities in the annulus K = <\t\ < 2}. Then 
the following bounds can be explicitly computed. 

1. the norm of the monodromy M corresponding to the circle {\t\ = 1} C JsT 
and its matrix logarithm B\ 

2. the pointwisc upper bound for norm of the fundamental solution X{t) with 
X{1) = E and its inverse X~^{t) in any smaller annulus K' , say, {| < \t\ < 

2 J ' 

3. the pointwise upper bound for ||W^(i)|| + ||W^~^(i)|| in the smaller annulus 
K'. 

In order to estimate the Laurent coefficients at all singular points after the 
surgery, it would be sufficient to find factorization (16.5) and supply the following 

bounds, 

1. maxtgDu ||i?o(OII + ll-f^o"^!^)!! (this would allow to estimate the norms of 
residues at all finite singularities). 
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2. the bound v for the order of the pole of 

3. maxfei,^ llC+,o(t)|l + \\C-.Q{t)\\ together with 

4. ^Yl^i^i l|C'±,i|| to majorize the norms of all Laurent coefficients at infinity. 

16.6. Quantitative matrix factorization. Unfortunately for our purpos- 
es, the known methods of constructing the Birkhoff-Grothendieck factorization 

(16.4), based on index theory for integral operators, do not allow for quantitative 
conclusions. Moreover, in some sense the problem admits no solution. The reason 
for this is the known instability of the partial indices z/i, . . . , an arbitrarily small 
variation (in the uniform norm) of the fimction W can result in a jump of the partial 
indices. This is clearly incompatible with existence of any bounds that would be 
continuous in the C'^-norm. 

However, if we give up with the uniquely defined Birkhoff-Grothendieck de- 
composition, then one can satisfy all the above conditions. 

Theorem 14 (Novikov and Yakovenko [NYOla]). A nxn-matrix function W{t) 
holomorphic and holomorphically invertible in the annulus 

K = {{l + 2e)-^ <\t\<{l + 2e)}, £ > 0, (16.8) 

and hounded together with its inverse there, 

\\W{t)\\ + \\W-'^{t)\\<R<+oo, tGK, (16.9) 

can be factorized as W{t) = HQ^H^{t) with the matrix functions Hq, satisfying 
the following conditions. 

1. iJo(i) is holomorphic invertible in the disk D'q = {\t\ < (1 + e)} and satisfies 
the inequality ||i?o(OII + 1 1 ^o"^ (01 1 ^ ^^^^ disk. 

2. Hoo{t) is holomorphic and holomorphically invertible in the complement 
D'^ = {(1 + < \t\ < +00} and both and have at most a 
pole of order v at t = 00; 

3. the coefficients C±^i of the Laurent expansions 

V 

P±{t)=Ht\t)=J2C±,if 

i=0 

are all bounded by R' in the sense of the matrix norm; 

4. the "regular parts" are bounded uniformly in D'^ so that \\H^ — P+{t)\\ + 
\\H--P.it)\\^R' there. 

The integer parameter v and real parameter R! can he expressed as explicit elemen- 
tary functions ofn,R ande, the "width" of the annulus K. 

16.7. Conclusion. Theorem 14 provides the last tool necessary to run the 

inductive proof and construct an explicit primitive recursive upper bound for the 
number of isolated roots of solutions to Fuchsian systems. The bounds can be 
explicitly written down in the sense that all primitive recursions describing them, 
can be extracted from the constructions above. However, this has not been done, 
among other things, because of the very excessive bounds that appear on this way. 
Besides nested inductive constructions, this occurs because of very excessive bounds 
on the lengths of ascending chains of polynomial ideals (sec §11). However, some 
very recent works of A. Grigoriev suggest that this crucial step can be considerably 
improved and instead of a tower of four stories, something like a double exponential 
bound can be achieved, at least for linear systems with rational coefficients. This 
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improvement may affect all other parts of the global construction as well, resulting 
in a considerably more realistic quasialgebraicity-type statements. 
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